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Chapitrel 

Variables, domes statistiques, 
tableaux, effectifs 


1.1 D definitions fondamentalesi.i.i La science statistique- M'ethodescientific 


- Etymologiquement: science de r'etat.- La statistique s'applique 'a la plupart des disciplines : agronorr 

biologiedemographissconomie, 

sociologie,linguistique,psychologie,... 1.1.2 Mesure et variable- On s’int'eresse 'a des unit'es statii 


des en- 

treprises, des rrTenages. En sciences humaines,on s’int’eresse dans la plupart des cas 'a un nombrefi 
d’unit’es.- Sur ces unities, on mesure un caract'ere ou une variable , le chiffre d’affaires de I’entreprise, 

le revenu 

du m’enage, l"'age de la p ersonne, la cat'egorie so cio-professionnelle d’une personne. On suppose c 
variable prend toujours une seule valeur sur chaque unite. Les variables sont d'esign'ses par simplic 
par une lettre ( X.Y.Z).- Les valeurs possibles de la variable, sont appellees modaiit'es - L’ensemble des val 
fait de pouvoir ou non ordonner les modaiit’es est parfois discutable. Par exemple :dans les categories soc 
Une variable est dite continue, si I’ensemble des valeurs possibles estcontinu. Remarque 1.1 Ces d'efinitic 


Exemple 1.1 Les modalels'de la variable sexe sont mascu//ne(tst)®17em/n/n (cosF). Le domaine 
de la variable est {M , F }. 


Exemple 1.2 Les modafc' de la variable nombre d'enfants par famille sont 0,l,2,i;'^S;.une va- 
riable quantitative diete'. 


1.1.4 S'erie statistique 

On app elle s'erie statistique la suite des valeurs prises par une variable 2fesutf'bfeBrwfcTtion. 
Le nombre d’unit'es d'observation esshnot' 

Les valeurs de la variable* sont eefc' 


Xi, ...,Xn. 


Exemple 1.3 On s'int'eresse 'ala variable ''etat-eB^railitalaserie statistiquedes valeurs prises par 
X sur 20 personnes. La codification est C: c'elibataire,M : mari'e(e),V : veuf(ve),D : divorc'ee.Le domaine d( 


MMDCCMCCC MC MVMVDCCC M lei, n= 20, xi =M,x 2 =M,x 3 =D,X4 =C,xs =C ,X20 =M. 1.2 Var 

I'effectif de la modalit'e xy. La frequence d'une moda li t'e est I 'e ff ectif divis'e par lenombre d’unit'es d'obs 


En langageR 

>X=c('Mari'e(e)', 'Mari 'e(e)7Divorc'e(e)7C'elibataire','C'elibataire', 'Mari 'e(e)','C'elibataire', 
'C 'elibataire', 'C 'elibataire', 'Mari 'e(e)','C 'elibataire', 'Mari 'e(e)','Veuf(ve)', 'Mari 'e(e)’, 
’Veuf(ve)','Divorc'e(e)','C 'elibataire', 'C'elibataire','C'elibataire','Mari'e(e)’) 
>Tl=table(X)> Vl=c(Tl)> data.frame(Eff=Vl,Freq=Vl/sum(Vl)) 


Eff FreqC'elibataire 9 0.45Divorc'e(e) 2 O.lOMari'e(e) 7 0.35Veuf(ve) 2 0.101.2.2 Diagrammeensed 


Le tableau statistique peut "etre repr'esent'e par un diagramme en barres ou en secteurs (ou camemb < 
ou piecharte n anglais) (voirFigures 1.1 et 1.2). Ceiibataire 



Fig. 1.1 - Diagramme en secteursEn langage R> pie(Tl,radius=1.0)En langage R> barplot(Tl)1.3 Va 


Celibataire Divorce(e) 


Marie(e) 


Veuf(ve) 


Fig.l.2-Diagramme enbarres 

La notation xi <: xi se lit xi pr'ec' ede X 2 .Si la variable est ordinale, on peut calculer les effectifs cumuT 

es: 

Nj =i 

k=ink. 

On a A/i = ni et Nj = n. On peut 'egalement calculer les fr'equen^sscunnur 

Fj = Njn = jk= 1 fk Exemple 1.5 On interroge 50 personnes sur leur dernier dipl'orne obtenu (variable Y ). 


La co dificationa 

'et'efaite selon le Tableau 1.1. On a obtenu la s'erie statistique pr'esent'ee dans le tableau 1.2. Finalemen 
obtientletableaustatistiquecompletpr'esent'edansleTableaul.3.Tab.l.l-CodificationdelavariableYT>ernierdi|: 
UEnlangageR>YY=c("Sd","Sd","Sd","Sd","P","P","P","P","P","P","P","P","P","P","P","Se","Se","Se","Se","Se",' 



Tab. 1.3- Tableau statistique complet 


x i 

nj 

Nj 

fj 

Fj 

Sd 

4 

4 

0.08 

0.08 

P 

11 

15 

0.22 

0.30 

Se 

14 

29 

0.28 

0.58 

Su 

9 

38 

0.18 

0.76 

U 

12 

50 

0.24 

1.00 


50 


1.00 



T2=table(YF)V2=c(T2)> data.frame(Eff=V2,EffCum=cumsum(V2),Freq=V2/sum(V2),FreqCum=cumsum(V2/ 
Eff EffCum Freq FreqCumSd 4 4 0.08 0.08P 11 15 0.22 0.30Se 14 29 0.28 0.58Su 9 38 0.18 0.76U 12 50 0.2- 


Figure 1.3). Sd 



Fig. 1.3- Diagrammeen secteursdesfr'equencesEn langage R> pie(T2,radius=l)9 


1.3.3 Diagramme enbarres deseffectifs 

Les effectifsd'une variablequalitative sont lasEtBriisau moyend'un diagramme en barres (voir Fi- 
gure 1.4). sd 



P Se Su u 

Fig. 1.4- Diagrammeenbarres deseffectifs 

En langage r> barplot(T2)1.3.4 Diagramme en barres des effectifs cumuTesLeseffectifs 


Figure 1.5). Sd 


8 



P Se Su U 


Fig. 1.5-Diagramme en barresdes effectifs cumul'eslO 



En langageR 

>T3=cumsum(T2) 

>barplot(T3) 


1.4 Variable quantitative disbe' 

1.4.1 Le tableaustatistique 

Une variable discr'etea un doma©rn©iahbrable. 

Exemple 1.6 Un quartier est compos'e des50gBS, et la variableZ repente le nombre de personnes 
par m'enage.Les valeursde lavariablesontl 1 1 11 


2 2 2 2 2 

2 22 233 3 3 3 3 

3 33 33 33 3 3 4 

4 4 4 4 4 444 4 5 

5 5 5 5 5 66 68 8 


Comme p our les variables qualitatives ordinales, on peut calculer les effestelKectifs cumeCles 
frequences, les frequences cumuTees. 'A nouveau, on peut construire le tableau statistique: 


xjnjNjfjFjl 5 5 0.10 0.102 9 14 0.18 0.283 15 29 0.30 0.584 10 39 0.20 0.785 6 45 0.12 0.906 3 48 0.06 0 



Fig. 1.6 - Diagramme en bfatonnets des effectifs pour une variable quantitative discr' 


1.4.2 Diagramme en b^atonnets des effectifsQuand la variableest discr'ete, les effectifs sor 


En langage R> plot(T4,type="h",xlab="",ylab="",main = "",frame=0,lwd=3)1.4.3 Fonction de repi 


esenie 

enFigurel.7,estd'efiniedeR dans[0,l]etvaut:F(x)=DDD 0 x<xiFj xj <x<x/+ ll xj <x.En langage R> plot(ecdf( 
discr'etes. Cependant, il est souvent int'eressant de proc'eder 'a des regroupements en classes pourfairec 


Fig. 1.7 - Fonction de repartition d'une variable quantitatiefeediscr' 


Exemple 1.7 On mesure la taille de 50 'el'eves d'une classe :152 152 152 153 153154 154 154 155 1551! 
50 0.20 1.0050 l.OOLe tableau re group e en classe est souv ent appel'e distribution group 'ee. On note, de 



-N j I'effectif cumeide la classe j, 

- fj la fr'equencede la classe j, 

- Fj la frequence cunedde la classe j. 


Enlangage R> S=c(152, 152, 152, 153, 153, 154,154, 154, 155,155, 156, 156, 156, 156, 156, 

+ 157,157,157,158,158,159,159,160,160,160,161,160,160,161,162, 

+ 162,162,163,164,164,164,164,165,166,167,168,168,168,169,169, 

+ 170, 171,171, 17 1,171)> T5=table(cut(S, breaks=c(151, 155, 159, 163, 167,171))) 

>T5c=c(T5)> data.frame(Eff=T5c,EffCum=cumsum(T5c),Freq=T5c/sum(T5c),FreqCum=cumsum(T5c/sum(' 
Eff EffCum Freq FreqCum(151,155] 10 10 0.20 0.20(155,159] 12 22 0.24 0.44(159,163] 11 33 0.22 0.66(162 


hauteur) repr'esente I'effectif. La hauteur hj du rectangle corresp ondant 'a la classe y'fflstpdonc donn' 
hj = n jc+j - c-j .o 



151.5 155.5 159.5 163.5 167.5 171.5 


Fig.l.8-HistogrammedeseffectifsEn langage R> hist(S,breaks=c(151.5,155.5,159.5,163.5,167.5,171.5),xla 


Si les deux denaiies classes sontagpfes,comme dans la Figure 1.9, la surface du dernier rectangle est 
'egale'ala surface des deuxetasirectangles de I'histogramme de la Figure 1.8. 



151.5 155.5 159.5 163.5 171.5 


Fig. 1.9- Flistogramme des effectifs avec les deu>ectemilasses ae^es 


En langage R> hist(S,breaks=c(151.5,155.5,159.5,163.5,171.5),xlab="",ylab="",main = "",xaxt = "n"; 




En langageR 

>y=c(0,0,cumsum(T5c/sum(T5c)),l) 

>x=c(148, 151.5, 155.5,159.5,163.5, 167.5,171.5, 175) 

> plot(x,y,type="b",xlab="",ylab="",xaxt ="n") 

> axis(l, c(151. 5, 155. 5, 159. 5, 163. 5, 167. 5, 171. 5)) 


16 



Chapitre2 

Statistique descriptive unwari' 

2.1 Param' etres de pOSition2.1.1 Le rnodeLe mode est la valeurdistincte correspond; 


eleve'; il est n(£b <m . 

Si on reprend la variable 'Etat civil’ , dontle tableau statistique est le suivant: 


X] nj fjC 9 0.45 M 7 0.35\/ 2 0.10D 2 0.1 0n= 20 lie mode est C: c'elibataire. Remarque 2.1- Le mode peut 


pondant 'a I’effectif le plus 'elev’e).2.1.2 La moyenneLa moyenne ne peut "etre d'efinie que sur une 


On peut aussi faire les calculs avec les valeurs distinctes et les effectifsef&nlEdiabldau: 


x j n i 

0 2 

1 3 

2 1 

3 1 

4 1 
8 


'x = 2x 0+ 3x 1+ 1 x2 +1 x 3+ 1x4 
8 

= 3+2+3 +4 
8 

= 1.5. Remarque 2.2 La moyenne n'est pas n 'ecessairementune valeurpossible. 


En langage RE=c(0,0,l,l,l,2,3,4)n=length(E)xb=sum(E)/nxbxb=mean(E)xb2.1.3 Remarques Slir 
X31 + X32 (/ = 3)18 



5. On peut exclure une valeur de I'indice. 


X; =x i +x 2 +x 4 +x 5. 


Propri'et'e 2.11. Somme d'une constante n/=i a= a+a+--+a 


=na (a constante). 


Exemple s/=i 3 = 3+ 3+ 3+ 3+ 3 = 5x 3 = 15. 


2. Mise en 'evidence ni=iaxi =an/=ix; (a constante). 


Exemple 3/=i 2x /= 2(1 +2 + 3) =2 x6 = 12.3. Somme des n premiers entiersn/=i i= 1+ 2+ 3+ ■■■+ n= n 



2.1.4 Moyenne gometrique 

Si xi > 0, on app elle moyeneensptrique la quanbt' 

1/n 1 In 1 n 

G= Xi =(x 1 XX 2 x ■■ ■ XXn) = exp- logXJ. 

La moyenne g'eom'etriques’utilise, par exemple, quand on veutcalculer la moyennesrifc.taux d'int ' 

Exemple 2.3 Supposons que les tau>ei£lirjl6ur 4 arses corescutives soient respectivement de 5, 10, 15, 
et 10%. Que va-t-onobtenir apr4 ans si je place 100 francs? 

- Apr'es 1 an on a, lOOx 1.05 = 105 Fr. 

- Apr'es 2 ans on a, 100 x 1.05 x 1.1 = 115.5 Fr. 

- Apr'es 3 ans on a, 100 x 1.05 x 1.1 x 1.15 = 132 . 825 Fr. 

- Apr'es 4 ans on a, 100 x 1.05 x 1.1 x 1.15 x 1.1 = 146 . 1075 Fr. 

Si on calcule la moyenne arithm'etiquedes faux on obtient 

'x = 1. 05 + l^LO -fr 1^ 5 + 1.10 


Si on calcule la moyenne g'eorrTetrique des taux, on obtientG = (1. 05 x 1. 10 x 1. 15 x 1. 10)i/4 

= 1. 099431377. 

Le bon taux moyen est bien G et non 'x, car si on applique 4 fois le taux moyen G aux 100 francs, on obtier 
100 Fr x G4 = 100 x 1.0994313774 = 146.1075 Fr. 


2.1.5 Moyenne harmoniquesi x; > 0, on appelle moyenne harmonique la quantit'eH = nn/=il/x;.ll e 
20b50= 20.8333/7,sa vitesse moyenne est done Moy= 40020.8333 = 19.2 km/h.- Si on calcule la moyenne 


Remarque 2.3 II est possible de montrer que la moyenne harmonique est teb^unesourfgale 'a la 
moyenneg'eom'etrique qui est toujeaesuiiref'oiffi'gale 'a la moyenne aratttropcie 


2.1.6 Moyenne pond 'er'eeDanscertains cas, on n’accorde par le rrT 

eme poids 'a toutes les observations. Par exemple, si on calcule 
la moyenne des notes pour un programerbecbgpn p eut pointer les notes daStudiant par le nombre 
decr'editsou parle nombred'heures dechaquecours.^idV / = 1, .. n sont les poids aesod'chaque 
observation, alors la moyenne pond'er'ee paefimiesfbad": 

~ xw i mm 

1 1 Wi 


Exemple 2.5 Supp osons que les notes soient pond'er'ees par le nomfaitdtsletEque les notes de 
T'etudiant soient les suivantes : Note 5 4 3 6 5Cr'edits 6 3 4 3 4La moyenne pond'er'ee des notes par les i 




2.1.7 La m "edianeLa m'ediane, not'ee xi/ 2 , est une valeur centrale de la s'erie statistique obtenue d 

- On trie la s'erie statistique par ordre croissant des valeurs observ'ees. Avec la s'erie observ'ee :3 2 1 
'epartition pour la valeur 1/2 :xi /2 =F-i(0.5).En langage R 21 



Fig. 2.1- M'ediane quand n est impair 



x=c(0 , 0 , 1 , 1 , 2 ,2 , 3)median(x)plot(ecdf(x),xlab="",ylab="",main = "",frame=FALSE,yaxt = "n") 


axis(2, c(0.0 f 0.25,0.50,0.75,1.00))arrows(-l,0.5,l,0.50,length=0.14,col="blue") 
arrows(l,0.50,l,0,length=0.14,col="blue")- Si n estpair,deuxvaleurssetrouventau milieudela s' 

erie (ici avec n= 8) 

0 01 1 2 234 

t t La m'ediane est alors la moyenne de ces deux valeurs :xi /2 = 1+22 = 1.5. La Figure 2.2 montre la 


ediane peut toujours 

"etre d'efinie comme I'inverse de la fonction de repartition p our la valeur 1/2 :xi /2 = F - i (0.5). Cep ei 


exactement'aun 'palier'. Fig. 2.2- M 'edianequand n estpair-i 



Enlangage Rx=c(0 ,0,1,1, 2, 2, 3, 4)median(x)plot(ecdf(x),xlab="",ylab="",main = "",frame=FALSE,ya 


En general on note 


*(D X(i) X(n) 

la s'erie ordoraipar ordre croissant. On app elle cffittiteca'don®e'la statistique d'ordre. Cette notation, 
tr'es usuelle en statistique, permotfchirda nadiane de maarrie tes syntbfique. 

- Si n estimpair 

X1/2 =X 

- Si n estpairx 

1/2 = 2 X (i) +x ( n 2+i) ■ 

Remarque 2.4 La m'ediane peuRetre ealsuKdes variables quantitatives et sur des variablesqualita- 
tives ordinales.2.1.8 QuantilesLa notion de quantile d'ordre p (o'u 0< p< 1) g'en' 


eralise la arfiane. Formellement un quantile esfedonn' 
par I’inverse de la fonction de Repartition : x P =F-i(p). 

Si la fonction de Repartition 'etait continue et strictement crofesfciftettiondu quantile serait sans 
'equivoque. La fonction de Repartition est cep endant disconfipBiiepEtiier".Quand la fonction de 
r'epartitionest par palier, il existe au moins 9 mani'eres diff'ereBfirardesdquantilesselon que I’on 
fasse ou non une interp olation de la fonction de Repartition. Nous pr'esentonsetrtiraafesjeaisrir 
ne faut pas s''etonnerde voir les valeurs des quantiles diff'erer I'eg'erement d’un logiciel statistique 'al'au 
- Si np est un nombre entier, alors x P = 12 x<np)+ x<np+i) .- Si rip n'est pas un nombre entier, alors x P = x< np 
statistique 12, 13, 15, 16, 18, 19, 22, 24, 25, 27, 28, 34 contenant 12 observations^ = 12).- Le premier qu: 


- Le troisieme quartile : Comme np = 0.75 x 12 = 9 est un nombre entier, ona 
_ *(9) +*(io) _ 25 +21 _ 


En langage Rx=c(12, 13,15, 16, 18, 19, 22, 24, 25, 27, 28, 34) 

quantile(x,type=2)Exemple 2.7 Soitlas'eriestatistique 12,13, 15, 16,18, 19, 22, 24,25, 27 contenantlOobs 


- Le premier quartile : Comme np = 0.25x 10 = 2.5 n'est pas un nombre entier, ona 

XI /4 = X( 2 .5) =X (3) = 15. 

- La m'ediane : Comme np = 0.5x 10 =5 est un nombre entier, ona 

xi /2 = 12 x<5) + x<6) = (18 + 19)/2 = 18.5. 

- Le troisi'eme quartile : Comme np = 0.75 x 10 = 7 .5 n’est pas un nombre entier, ona 

X3/4 = x< 7.5 ) = X(8) = 24. En langage Rx=c(12,13,15,16,18,19,22,24,25,27)quantile(x,type=2)2.2 Para 


2.2.3 La variance 

La variance estla somme des earcdesecarts 'a la moyennectieipar le nombred’observations: 
s 2 x = ^ (Xi - 

Th'eor eme 2.1 LavariancepeutaussiBi're 



(2.1) 


D'emonstration S 2 x = In n/=i(x/ - "x )2 =1 


(xf -2xi~x + 2 )x 

l=i 

= In ni=l X2i -2 In ni=l Xi'/Xn+l n n „ i n 

- '1= - xf - 2"* X; + "3e 

;=i n M n ;=i 

= In ni= i X2/ - 2 x x + X2 — In m=i X2\ - "x 


La variance peut 'egalement'etre d 'efinie ' a partir des effectifs et des valeurs distinctes: 

S2x = lnjj=inj(xj - "x) 2 .La variance peut aussi s''ecrire S2x = lnjj=injX 2 j - "x 2 .Quand on veut estimer une va 


s'electionn'ee au hasard) de taille n, on utilise la variance "corrig'ee" divis'ee par n-l.S 2 x = ln-1 n/=i(x/- 
2 + 3+4+4+5 + 6+7+98 = 5,25 



s*= - (X; - 

11 i= 1 

= 1 8 (2 - 5f + (3 - 3) + (4 - 3)+ (4 - $1+ (5 - $+ (6 - 5*)+ (7 - 3^+ (9 - 3) 
= 18 [9 + 4 + 1 + 1 + 0 + 1 + 4 + 16] 

= 368= 4. 5. On peut'egalement utiliser la formule (2.1) de la variance, 


la moyenne n'est pas un nombreentier). 

S2x = In ni = 1 X2i - ~X 


ce quinecessitemoins de caldsUrtout quand 


= 18(4 +9 + 16 + 16 + 25 + 36 + 49 + 81) - 25 


= 2368 -25= 29. 5-25=4.5. EnlangageR> x=c(2,3,4,4,5,6,7,9)> n = length(x)> s2=sum((x-mean(x))''2)/n 



2.2.5 L''ecart moyen absolu 

L ' ' ecart moyen absolu est la somme des valeurs abscriaestslea'la moyenne (iie'ipar le nombre 
d'observations: 

Qn oy = \Xi ~ X \. 

i= 1 

2.2.6 L'ecartm'edian absolu 

L' 'ecart m 'edian absolu est la somme des valeurs absBdaetedeslaediane divee par le nombre 
d'observations :e 

_ 1 n 

med ~ Xj —X x /2 ■ 

n i = 1 

2.3 MomentSD' efinition 2.2 On appelle moment 'a I'origine d’ordre rGN le param' 

etre 


m r = Inn 

- ;= ix ■. 


D 'efinition 2.3 On appel le moment centr'e d’ordre r E I Miasparam' 

mr = In ni= i (xi - ~x)r . Les moments g'en'eralisent la plupart des param' etres. On a en particulier- mi = ~x 


etrie 

et rapiatissement.2 .4 Param' etres de forme2.4.l Coefficient d'asym 'etrie de F 

normalis'e par la distance interquartile :Ay = X3/4 +xi/4 -2xi/2X3/4 -xi/4 .27 



2.4.3 Coefficient d'asyokrie de Pearson 


Le co efficient cTas§trie de Pearson est beisur une comparaison de la moyenne et du rmfe det 
standardise paadcirt-typ e: 



Tous les co efficients d'as^tnre ontles rremes propeites, ils sont nuls sila distribution est syetrique, 
n'egatifs si la distribution est ailong gauche (left asymmetry), et positifs si la distribution ®si ailong' 
droite (rightasymmetry) comme medfetns la Figure 2.3. 



Fig. 2.3-AsyrrTetrie d'une distribution 


Remarque 2.6 Certainesvariablessonttoujours tr'esasyrrTetriques'adroite, comme les revenus,les tallies 
des entreprises, ou des communes. Une m'etho de simple pour rendre uneatoiqaitelEsisrEiiis'te alors 
'a prendre le logarithme de cette variable. 2. 5 Param'etre d'aplatissement (kurtOSi 
courtes.Dans la Figure 2.4, on pr'esente un exemple de deux distributions de rn'eme moyenne et de m'em 



Fig. 2.4- Distributions m'esokurtiqueet leptokurtique28 



2.6 Changementd'origine et d'ueit 

D efinition 2.4 On appel le changement d'origenstltjrpc'onsistant 'a ajouter (ou soustrairejetaien' 
quantit'e aE R 'a toutes les observations 

y-, = a+Xj, i = 1, ...,n 

D efinition 2.5 On appel le changementlCbpefb'tionconsistant 'a multiplier (ou diviser) pasrrtBem'' 
quantit'e bE R 'a toutes les observations 

y =bxi, i = 1, n. 

D efinition 2.6 On appel le changedlerigine etd'unibVoperation consistanta multiplier toutes les 
observations par la m'eme quantit'e bE R puis 'a ajeoteiqliaamntitaE R 'a toutes les observations: 

yi = a+ bn, i = 1, n. 

Th'eor eme 2.2 Si oneffectueun changementd'origine ebdfurrih' variable x, alors sa moyenne est 
affect 'ee du m'eme changement d’origine et d’unit'e. D'emonstration Si y/ = a+ bxi, alors’y = In n/=i (a 


i=lX 


= a+ b x. 


* 

Th eor eme 2.3 Si on effectueun changement d’origine et d’unit'esur une variable x, alors sa varianceest 
affect 'ee par le carr'e du changement d'unit'eet pas parle changement d'origine. 

D'emonstration Si y; = a+ bxi, alorss 2 y = In n/=i (y< - ~y ) 2 = In ni=i (a + bxi -a - b ~x)i = bi In ni=i (xi - 


Remarque 2.71. Les param'etres de position sont tous affect'es par un changement d'origine et d'unit'e. 


d'origine. 3. Les param'etres de forme et d’aplatissement ne sont affect'es ni par un changement d'unit'e r 
, xha+i, xnA+2, .. ., Xn-i, Xnobservations deGs .On d'efinit les moyennes des deux groupes : 29 



1 nA 

- la moyenne du premier group @= ■*- x-,, 

^ i=1 n 

- lamoyenne du deuecrhe group % — x,. 

115 i=n A+ 1 

La moyenne g'en'erale est une moyeneie^padla taille des groupes des moyennes des deux group es. 
En effet “x= 

- x, + x, = 1 (n A ~X +n B 'X )■ 

n 1=1 1=04+1 n 


On p eut 'egalement d'efinir les variances des deux group es: 

i n a 

- la variance du premier group e=sa- ( Xi - » 2 , 

nA /= i 

- la variance du deuxi' eme group e sib ="l _ , 2 

— (X/ - &■) . 


Th'eor'eme 2.4 (de Huygens)Lavariancetotefmab'par 

S2x = In ni= 1 (Xj - ’X) 


se d'ecompose de la mani'ere suivante :s2x = has2a + ne S2Bnvariance intra-groupes + m( xa - 'x)2 +n 

b Cx - 

nv aria ne e inter-group es 


D emonstration S 2 x = In ni=i(xi - 'x )2 = ln/w=i(x/ - 'x )2 + ni=n A +i(xi - 'x )2 (2.2) 


On note quen4/=l(x/ - x)2= nM=l(Xi - XA + X4 - x)2= 04/=l(X/ - X4)2 + r)Ai=l{~ XA - 'x)2 +2 04/= 1 (x/ - -XA)('X> 



2.8 Diagramme en tigesetfeuilles 

Lediagrammeentigeset feuillesou Stem and learfiagrame stune manare rapide de pasenter une 
variable quantitative. Par exemple, si I'on aria st'atistique ordoraisuivante: 

15,15, 16, 17, 18, 20,21, 22, 23,23, 23, 24, 25, 25, 26, 

26,27, 28, 28, 29, 30,30, 32, 34,35, 36, 39, 40, 43, 44, 
la tige du diagramme sera les dizaines et les feuilles seron4tei>iiindtt5tient le graphique suivant. 

The decimal point is 1 digit(s) to the right of the| 

1 | 556782 | 0123334556678893 | 00245694 | 034Ce diagrammepermet d’avoirune vuesynth'etiquede lad 


Evidemmentles tiges peuvertifre 

d'efinies par les centaines, ou des millers, selon I'ordre de grandeur destaicki&iable' 

En langage R## Diagramme en tige et feuilles#X=c(15,15, 16, 17, 18, 20,21, 22, 23,23, 23,24,25, 25,26, 26,2 1 


qui p ermet de repr'esenter la distribution d’une variable. Ce diagramme est comp os'e de Un rectangle 

pondant 'a la rrTediane.- Ce rectangle est compl'et'e par deux segments de droites.- Pour les dessiner, on 
extr'emes”. Exemple 2.9 On utilise une base de donn'ees de communes suisses de 2003 fournie par I’Offic 
'etriedeladistribution,ausenso'uilyabeaucoupdepetitescommunesetpeudegrandescommunes.Legraphiquen 


1000 


2000 


3000 


4000 


5000 


Fig. 2.5 - Bcfites 'a moustaches p our la variable sup erficie en hectares (HApoly) des communes du cantor 
de Neuchfatel# 'Etape 1: installation du package sampling 


# dans lequel se trouve la basede donn'eesdes communes beiges 

# choisir "sampling"dansla liste 

utils:::menulnstallPkgs()# Etape 2: charge le package sampling# choisir "sampling"dansla I isteloca l( { pkg < 


+ if(nchar(pkg)) library(pkg, character.only=TRUE)}) 

# Utilisation des donn'eesdata(swissmunicipalities)attach(swissmunicipalities)# boxplot de la selection de 


Exemple 2.10 On utilise une base de donn'ees beiges fournie par I'lnstitut National (b elge) de Statistique 
contenant des infirmations sur la population et les revenus des p ersonnes physiques dans les communes. ( 
s’int'eresse 'a la variable "revenu moyen en euros par habitant en 2004" pour chaque commune (variable 
rageincome) et I’on aimerait comparer les 9 provinces b elges : Anvers, Brabant, Flandre occidentale, Fland 
orientale,Flainaut, Li'ege, Limboug, Luxembourg, Namur. La Figure 2.6 contientles bo"ites'amoustachesdec 
noms des provincesb=list("Anv."=averageincome[Province= = l],"Brab."=averageincome[Province= = 2],"FI. 



Brab. Fl.c 


Fl.or. Hainaut Liege Limb. Lux. Namur 


Fig. 2.6 - BcTites 'a moustaches du "revenu moyen des habitants" des communes selon les provinces belgi 


"Namur"=averageincome[Province==9])boxplot(b)ExerciceSExercice2.1 On p'ese les 50 'el'eves d’ 


esunes dans le tableau 

suivant : 43 43 43 47 4848 48 48 49 4949 50 50 51 5152 53 53 53 5454 56 56 56 5759 59 59 62 6263 63 (. 
classes. Exercice 2.3 33 



1. Montrez que 


2. Montrez que 

3. Montrez que, si x; >0, 




1 

2n2 


(X; 


Sx <E t 


n-1 
2 n 


1 n 



|x/ - x| 


■xj) 2 . 


2~x. 


34 



Chapitre3 

Statistique descriptive biorari 


3.1 S eriestatistique bivari ee 

On s'int'eresse 'a deux variables x et y. Ces deux variables ssmis nete'n unafe'd'observation. 
Pourchaque unit'e,on obtient done deux mesureErlfflsba'tistique est alors une suite de n couplesdes 
valeurs prises par les deux variables sur chaque individu :(xi, yi), (x;, yi), (x 


Chacune des deux variables peut "etre, soit quantitative, soit qualitative. On examine deux cas. 

- Les deux variablessont quantitatives.- Les deux variablessont qualitatives.3.2 Deux variables 


r'eels) peut toujours^etre repr'esent'e comme un point dans un plan(xi,yi),.... / (x/ / y/),....,(xn / yn).Exemple 3, 




155 160 165 170 175 180 185 190 

taille 

Fig. 3.1- Lenuagedepoints 


3.2.2 Analyse des variablesLes variables x et y peuvent'etre analys'ees s'epar'ement. On peut c 


moyennes et les variances : x = In n)=i x/, S 2 x = In n/=i(x/ - "x) 


y= In ni=i yi, S 2 y = lnn/=i(y/- y) 2 . Ces param'etres sont appel'es param'etres marginaux : variances ma 


ecarts-types 

marginaux, quantiles marginaux, etc. ..3. 2. 3 CovarianceLa covariance est d'efinie sxy = lnn/=i(x; - ~x 




Demonstration 


Sxy = ^ (Xi-'xMy-y) 

= n (XiYi -yi'x yxi + yj 



1 - Xiyi - > - >+ X 

/=i 

= In n/=i xy 

/ - >" 


* 


3.2.4 Corr'elationLecoeffic/ent de corr'elatione stlacovariancedivis'eeparlesdeux' 

ecart-typ es marginaux: 

rxy = sxysx Sy .Le coefficient de determinati on es t le carr'e du co efficient de corr'elation: 


r2xy = S2xyS2x S2 y .Remarque 3.2- Le coefficient de corr'elation mesure la d'ependance lin'eaire entre deu> 


cependant avoir une d'ependance non-lin'eaire avec un coefficient de corr'elation nul. 3.2.5 Droite de 
commet (voir Figure 3.3) en utilisant la droite de degression pour pr'edirey/ 'a partir de x/. Les r'esidus pei 



poids 



Fig. 3.2 - Exemples de nuages de points et co efficientsattaioarr' 



155 160 165 170 175 180 185 190 


taille 


Fig. 3.3 - Le nuage de points, le r'esidu38 


En langageR 

# Graphique avec le r'esidus 
plot(taille,poids) 

segments(158,a + b*158,190,a + b*190) 
segments(180,a + b*180,180,96,col = "red") 

#text(178,90,expression(e)) 

text( 178.7, 89.5,"i")#arrows( 180, a + b*180,156,a + b*180, col = "blue",length =0.14) 


arrows(180,60,180,a + b*180,col = "blue", length =0.14) 
arrows(180,96,156,96,col = " blue", length =0.14) 

#text( 154.8, 86, expression(y))text(155.5, 85.5," i")#text(154.8,97,expression(y))text( 155.5, 97.8,"*")text( 155 


es qui consiste 'a 

chercher la droite qui minimise la somme des carries des r'esidus :/W(a, b) = m=i er\ = n/=i (y/ -a - bxi)i .Th 


es sont donars' par: 

b = sxyS 2 x et a= ~y - b'x. D emonstration Le minimum M ( a,b ) en a,b s’obtient en annulant les d'eriv'ees 


et b. □□□□□□□□□ dM(a,b)da =- m = i 2 (y/ -a- bxi) = 0 dM(a,b)db =- m= i 2 (y; -a- bxi) x, = OOn obtient un sy: 


La premereequation montre que la droite passe par le pojf)t Onxqbtient 
a =y - b~x. 

En rempla.cant a igad b'x dans la secoeqlaation, ona 


1 

n 


XiYi -( y-b'x ) 



1 - - u 1 2 -2 

m *- b n,.*-* 

=S xy -b^ x 
= 0 , 


ce qui donnes 
Done b= 


xy -b£ = 0 . 

Sxy 
S 2 x ' 


On a done identifi'e les deux param'etresDDDDD b= sxyS2x (la pente)a = y-b x= y - sxyS2x 


x (la const: 


On devrait en outre v'erifierqu’il s'agit bien d'un minimum en montrant querliaeoct'secondes sont 
positives.* 

La droite de degression est done y=a+bx= ~y— sxyS2x ~x+ sxyS2xx,ce qui peut s''ecrire aussi y- ~y = sxys 2; 



Remarque 3.3 Ladroitede r'egressionde y en x n’estpasla irTernequela droiteder'egressionde x en y .40 



3.2.6 R esidus et valeurs ajeefc 

Les valeurs ajuebs sont obtenuesau moyende la droitecjEeBsion: 
y f =a + bx. 

Les valeurs ajust'ees sont tenons' desyealisesau moyen de la variablex et de la dro&^retesrdn 
de y en x. Remarque 3.4 La moyenne des valeursajust' 

ees esfegale 'a la moyennedes valeursafceagirvEn effet, 

In ni= 1 y* i n i n 

i= - (a + brf = a +b h x, = a+ b'x. 

Or, "y = a+ b'x, car le point fx, 'y) appartient 'a la (dgatesiateir' 

Les r'esidussont les diff'erencesentre les valeeesBafefeta' valeurs a^iestde la variabtepdendante. 
e \ = y, 

Les r'esidusrepr'esententla partie inexpliqu'eedes y/ parla degiitesdisn: 

Remarque 3.5- La moyenne des r'esidus est nulle. En effetlnn/=i e/ = lnn/=i(y; - y*i )= "y - ~y = 0. 


- De plus, ni = i x/e/ = O.La d demonstration est un peu plus difficile. 3. 2. 7 Sommes de carries et var 



D efinition 3.4 On appel le somme des ds&esidus (ou&siduel le) la quaetit ' 


SC res = ef. 

D efinition 3.5 Lavariancesiduel le est la variance dekliu's. 

4 _ _ 1 ' ^ 

n n j=1 

Note : II n'est pas n 'ecessaire de centrertesdussur leurs moyennes pour calculer la varianaeja 
moyenne des r'esidus est nulle. Th'eor'eme 3.3 SCror 

=SC regr +SC res- 

Demonstration SC tot = ni=i ( y; - ~y) 2= ni= 1 (y; - y*; + y*/ - ~y) 2= m= 1 (y/ - y*; )i + ni= 1 (y*/ - 'y)2 +2 


H (y -y*)(y*-y) 


= SC res + SC regr +2 n/=i(y/ - y*; )(y*/ - y).Le troisi'eme terme est nul. En effet,n/=i(y; - y*; )(y*; - y) 



Demonstration 



Th eor'eme 3.5 La variance r'esidueiie peut 'egalement s"ecrires 2 e = S2y(l- n),o'u n est le coefficient < 

variables qualitatives3.3.1 Donn ees observeessi les deux variables x et y sont qua 


chacune des deux variables prend comme valeurs deesnpjdStSt'ives. 
Lesvaleursdistinctesde x et y sontnois resp ectivement 


*1 >f * 



3.3.2 Tableaude contingence 

Les donn'ees observ'ees pedii^nigroigE'ssouslaforme d'un tableau de contingence 



yi 

... y k .. 

■ y K 

total 

xi n 

ii 

nik ■■ 

niK 

n i. 

xj nj 

i ■■■ 

njk 

njK 

nj. 

xj nj 

i ■ ■ 

■ njk 

■ n/K 


total 

n. i ■ 

■ ■ n.kn 

.K 

n 


Les nj. et n.k sont app el'es les effectifs marginaux. Dans ce tableau, 

- nj. repr'esente le nombre de fois que la mo dalit'e xfc.appara'' 

- n.k repr'esente le nombre de fois que la mo dalit'e yk appara'it,- njk repr'esente le nombre de fois qu€ 
On a les relations jj=injk = n.k, pour tout k= l,...,K,Kk=injk = nj., pour tout j= l,...,J,etjj=inj. =Kk=in.k =jj=iKk=ir 


Le Tableau 3.1 reprend le tableau de contingence.Tab. 3.1 - Tableau des effectifs nyfcBleu Vert Marron TotalH< 


f k = nk ,k = l K. 

n 

Le tableaudes frequences est 



yi ■■ 

■ y k ■■ 

■ y K 

total 

Xi 

fn ■■ 

■ fi k ■■ 

■ flK 

fi. 

Xj 

fji ■■ 

■ f Jk ■■ 

■ fjK 

h 

X] 

fji " 

■ fjk ■■ 

■ fjK 

fj- 

total 

f. i ■■ 

• f. k 

f.K 

1 


Exemple 3.3 LeTableau3.2reprend letableaudesdiquences. 


Tab. 3.2- Tableau des fequences 


Bleu Vert Marron 

Total 

Homme 0.050.25 0.10 

0.40 

Femme 0.100.30 0.20 

0.60 

Total 0.15 0.55 0.30 

1.00 


3.3.4 Profils lignes et profils colonnesuntableau de contingences'interpr'etetoujours en comp 

equences 

en colonnes (appel'es aussi profits lignes et profils colonnes ). Les profils lignes sont d'efinis parfm = njknj. = 



Tab. 3.4 - Tableaudes profils colonnes 



Bleu 

Vert 

Marron 

Total 

Homme 

Femme 

0.33 

0.67 

0.45 

0.55 

0.33 

0.67 

0.40 

0.60 

Total 

1.00 

1.00 

1.00 

1.00 


3.3.5 Effectifs theoriqueset khi-earr 

On cherche souventune interactionentre des lignes et des colarmbep entreles variableSour 
mettre en 'evidence ce lien, on construit un tableau d'effeantiic^fels quiepie'sente la situation o'u les 
variables ne sont pas li'ees (ind'ep endance). Ces asffkiptds febnt construits de la leransuivante: 


Les effectifs observ'es njk ont les irfemes marges que Ieseeffqctel3$th' 
Enfin, les 'ecarts 'a I’ind'ependance sont d'efinis pareyfc = rijk - n* 


jk. 


- La d'ep endance du tableau se mesure au moyen dedkfelri-iaBpta'r 

X2 obs = Kk=ljj=l ( njk - n*jk)2n*jkj^fgl3=lj 


(3.1) 


- Le khi-carr'e p eut "etre normalis'e p our ne plus d'ep endre du nombre cOmbdefiiviatierfsbi- 
deux par : (p 2 = x 2 °bsn .Le <pi ne d'epend plus du nombre d'observations. II est possible de montrer que<p 2 < 
tableau ait le rrfeme nombre de lignes quede colonnes). Exemple 3.5 Le Tableau 3.5 reprend le tableau de 


Tab. 3.6 - Tableau deecarts 'a I'iepl endanc#e 



Bleu 

Vert 

Marron 

Total 

Homme 

-2 

6 

-4 

0 

Femme 

2 

-6 

4 

0 

Total 

0 

0 

0 

0 


Tab. 3.7 - Tableaudes^ fc /n ; £ 


Bleu 

Vert 

Marron 

Total 

Homme 0.33 

0.82 

0.67 

1.82 

Femme 0.22 

0.55 

0.44 

1.21 

Total 0.56 

1.36 

1.11 

3.03 


- Le khi-carr'e observ'e vaut * 2 obs = 3.03.- Le phi-deux vaut qn = 0.01515.- Comme le tableau a deux li< 


au <j92 .- On a V = 0.01515. La d'ep endance entre les deux variables est tr' 

es faible. 


egal 


Enlangage Ryeux= c(rep("bleu",times=10),rep("vert",times=50),rep("marron",times=20), 


rep("bleu",times=20),rep("vert",times=60),rep("marron",times=40))sexe= c(rep("homme",times=80),rep(" 
Tableaudecontingence: effectifsnyfcNiveaud'instruction Statutprofessionneldufilsdufilsparrapport parrapporta 



Tab. 3.9 - Tableau des fequence^ 


X\Y 

Pluseleve' Egal inlerieur 

total 

pluseleve' 

egal 

inferieur 

0.322 0.231 0.147 

0.055 0.079 0.058 

0.017 0.038 0.053 

0. 700 
0. 192 
0. 108 

total 

0.394 0.349 0.257 

1. 000 


Tab. 3.10- Tableau des profils lignes 


XirPlu$' 


eleve' Egal 


inferieur| 


total 


plus 'elev 
'egal 0. 
inf'erieu ' 


'e 0.460 
288 
0. 156 


0. 330 0. 210 1 

0. 413 0. 300 1 

0. 356 0. 489 1 

0. 349 0. 257 | 1 


Tab. 3.11 - Tableaudes profils colonnes 

X\Y P\ui 'elev'e Egal inf' erieur total 
ifgus 'elev'e 0.817 0.662 0.570 0. 700 

'egal 0.140 0 .228 0.224 0.192 

inf'erieu - 0.043 0 .110 0 .206 0. 108 


Tab. 3.12 - Tableau des effectifs th'eoriques n*jkX\Y Plus 'elev'e Egal inf'erieur totalplus 'elev'e 114.72 
1.24 12.52 *2 obs = 25.810n a done *2 obs = 25.81p2 = x* obsn = 25.81416 = 0.062\/ = <p 2 min(/ - 1,/C- 1) = 



ExerciceSExercice 3.1 Laconsommationdecr' 

ernes glaaes par individus^temesuee pendant 30epio des. L'ob- 
jectif est d'etermin'e si la consommafjoamd'de la terwiatu reJ.es doniees sont dans le tableau 3.15. 
Onsaiten outreque Tab. 3.15-Consommationde cr' 


ernes gl®es 


consommation y temp'eratur 

ex consommatiortf£m|£faturex 

consommationytemferaturex 

386 41 286 

28 

319 

44 

374 56 298 26 


307 

40 

393 63 329 32 


284 

32 

425 68 318 40 


326 

27 

406 69 381 55 


309 

28 

344 65 381 63 


359 

33 

327 61 470 72 


376 

41 

288 47 443 72 416 



52 

269 32 386 67 437 



64 

256 24 342 60 548 



71 


ni=i yi = 10783, ni=i xi = 1473,n/=/ yn = 4001293, ni=i X 2 i = 80145,™=/ x/y/ = 553747,1. Donnez les moyenne: 


2. Donnez la droite de degression, avec comme variable d'ep endante la consommation de glaces et cor 
variableexplicativelatemp'erature.3. Donnezlavaleurajust'eeetler'esidupourlapremi'ereobservationdutable 
'etationdececoefficient.3. D 'eterminezladroiteder'egressionY enfonctiondeX.4. 'Etablissez,surbasedumod' 


Exercice 3.3 Conskfons unechantillonde 10 fonctionnaires (ayant entre 40 et 50 ans) d'unermnist' 

Soit X le nombre d’aees de service et Y le nombre de jours d'absence pour raison de maladie (au cours de 
I'ann'ee peltoite) dtermie'pour chaque p ersonne appartenanteckaMiflon. 


Xi 

2 

14 

16 

8 

13 

20 

24 

7 

5 

11 

JL 

3 

13 

17 

12 

10 

8 

20 

7 _ 


8 


1. Repr'esentez le nuage de points. 

2. Calculez le co efficient deedathdn entre X et Y. 

3. D'eterminezr'equationdela droiteei^cetssionde Y en fonction de X. 

4. D'eterminez la qualit'e de cet ajustement. 

5. 'Etablissez, sur base de ce mo d'ele, le nombre de jours d'absence p our un fonctionnaire ayant 22 an 
service. 50 



Chapitre4 


Th'eorie des indices, mesures 

d'inegalit04.i Nombres indices4.2 Definitionun indice est ia 


efe'rence. Prenons I'exemple du tableau 

4.1 contenant le prix (fictif) d’un bien de consommation de 2000 'a 2006. Le temps varie de 0, 1, 2, .. ., 6e 
est consid'er'e comme le temps de r'ef'erence par rapp ort auquel I'eidice est calcul' 


Tab. 4.1- Tableau du prix d'un bien de co nsommation de 2000 'a 2006 
ann'ee t prix pt2000 0 2.002001 1 2.302002 2 2.402003 3 2.802004 4 3.002005 5 3.502006 6 4.00L'indice 
142.864 66.67 76.67 80.00 93.33 100.00 116.67 133.335 57.14 65.71 68.57 80.00 85.71 100.00 114.296 5 


4.2.1 Proprietes des indices 

Consioferons un indice quelconque / (t/0) . On dit que cet indeiiepl®®P'i'opaies de 
-r'eversibil&si I (t/0) = l(ft)x 

- identit'e si /(t/t) = 100, 

- circularit'efou transitiefcsi I ( t/u ) x I ( u/v ) = 100 x / (t/v). 

II est facile de montrerque ces quatre petqsrBont satisfaites pour un indice simple. 

4.2.2 Indices syntfet'iques 

Quand on veutcalculer un indice 'apartir de plusieurs prig, problame devient sensiblement plus 
compliqu'e. Un indice synth'etique est une grandeur d’un ensemble de biens par rapp oide' dome ann' 
r'ef'erence. On ne peut pas construire un indie&qagoTih additionnant simplement des indices simples. 
Ilfaut,en effet,tenir compte des quaielatach^es. 

Pour calculer un indice de prix de n biens de consomatiqtietes de 1,2, .. ,,n, on utilise la notation 
suivante:- pti repr'esenteleprixdu biendeconsommation / autemps t, 

- qn repr'esente la quantit'e de biens / co®eraramtrernps t. 

Consid'erons par exemple le Tableau 4.3 qui contient 3 biens de consommation et pour lesqitels ont conna' 
les prix et les quantities achet'ees.Tab. 4.3 - Exemple : prix et quantities de trois bien pendant 3 ans 


Temps 

0 12 

Prix ( poi ) Quantit'es ( qi 

>/ ) Prix ( pn ) Quantises ( 

q i i ) Pri #Xp Quantifes (g,) 

Bien 1 

100 14 150 10 200 


8 

Bien 2 

60 10 50 12 40 


14 

Bien 3 

160 4 140 5 140 


5 


II existe deux rrTetho des fondamentales p our calculer les indices de prix, I’indice de Paasche et I'indice 
Laspeyres.4.2.3 Indice de LaspeyreSL'indice de Laspeyres, est d'etini parL(t/0) = 100 x n ;= i qoiptini- 
woili(t)ni=i woi = ni= i po/Qo/100 x p H po/n/=i poiqoi = 100 x ni= i qoiptini=i poiqoi .L’indice de Laspeyres ne poss'ede 



Exemple 4.1 Sionutiliselesdonrees du tableau 4.3, les indices de Lasp eyres sont les suivants 


L(l/0) = lOOx- 
L( 2/0) = 100 x- 
L(2/l) = lOOx 


qb/pi; . 
'1=1 QaiPoi 


x — 1S0 + 10X50+4X 14Q 
14 x 100 +10 x 60+ 4x 160 ' 

14x 200 + 10 x 40+4 x 14() 


/=i QaiPoi 
l OiiPzi . 


14 x 100 +10 x 60+ 4x 
lOx 200 + 12 x 40+5 x 


10 x 150 +12 x 50+ 5x 140 


4.2.4 Indice de PaascheL'indice de Paasche, estd'efini parP(t/0) = lOOx 


i=l QtiPti 
i=i QiPoi' 

On utilise, pour le calculer, les quantities qa du temps par rapp ort auquel on veut calculer I'indice. 

L’indice de Paasche peut aussi "etre pr'esent'e comme une moyenne harmoarrapjdqsdjroiJces 
simples. Soient I'indice simple du bien / : //( C/0 ) = 100 x pt/po; ,et le poids wti correspondant 'a la recette tot 


etes par 

les recettes au temps t :P (t/0) = m= i wtini=i wti/li(t/0) = m= i pnqtini= i ptiqti poiiooxpti = 100 x n/=i qtiptini= i qtipoi .1 


difficile 'a calculer que I'indice de Laspeyres, car on doit conna^itre les quantities pour chaque valeur de t. 
FisherL'indice de Laspeyres est en g^n'eral plus grand que I’indice de Paasche, ce qui peut s'expliquer p 



Fisher a prop©sf utiliser un compromis entre I’indice de Paasche et de Lasp eyres en calculant simpleme 
la moyenneeaphetriquede ces deux indices 

F(t/ 0)= L{t/ 0) x P [t/ 0). 

L’avantagedel'indice de Fisher estqu'il jouit de la pelepdie reversibilS.' 

Exemple 4.3 Sionutilisetoujourslesdonaes du tableau 4.3, les indices de Fisher sont les suivants: 

Fi 1/0)= L(l/0) x P (1/0) = 115.3242, 

F(2/0)= L(2/0) x P (2/0) = 129.2052, 

F(2/l)= L(2/l) x P (2/1) = 111.7715. 


4.2.6 L'indice de SidgwickL'indice de Sidgwick est la moyenne arithm'etique des indices de Paasch 
S (t/0) = Ljt/ 0) 4 - P jt/ 0 ) 


4.2.7 Indices cha'inesLe d'efaut principal des indices de Laspeyres, de Paasche, de Fisher et de Side 

edent 

pas la propri'et'e de circularit'e. Un indice qui poss'ede cette pro|arn#tffi(datieppoeJr construire 
un indice cha^ine, avec l'indice de Lasp eyres, on peut faire un produit d'indice de Lasp eyre annuels. 

CL (t/0) = 100 x Ljt/t - 1 )1 00 x Ljt - 1/ t - 2)1 00 x- x ^1)100 x L(l/0) 

Pour calculer un tel indice, on doit 'evidemment conna'itre les quantises pour chaqueLtokiiiaede t. 
suisse des prix 'a la consommation est un indice cha'ine de Laspeyres. Exemple 4.4 En utilisant encore les 


suivants : CL(l/0)= L(l/0)= 119.6970,0.(2/1)= L(2/l)= 113.5714,CL(2/0)= L(2/l)xL(l/0)100 = 135.9416.4 


4.3.2 Courbe de Lorenz 


Plusieurs indices dbgalife'sont les 'a la courb e de Lorenz. On note 


Xi , ,x ... ,x 

les revenus des n individus de la populatiaxitiee. On notegalement 


*<« ft *». 

la statistique d'ordre, c'est-'a-diiffirfexle revenusisipar ordre croissant. 

Notons maintenant qi la proportion de revenus par rapp ort au revenu total qeresif gmjgviidus 
ayant les plus bas revenus,ce qerinif' 




" i i= 1 *(/> 

avecc * =0 etc ?" 


= 1 . 


La courb e de Lorenz est la representation graphique de la fonction qui 'a la part des individus les moin 
riches asso cie la part y du revenu total qu'ils p er.cBfofBpkecismentJa courb e de Lorenz relie les 
points ( i/n , qi) pour / = 1, .. ., n. En abscisse, on a done une proportion d'individuesc(tBEH5ordre de 
revenu, et en ordonn'ee la prop ortion du revenu total re.cu par ces individus. 


Exemple 4.5 On utilise une enqu'etern'enage sur le revenudagiis>DrcteE'Philippines apqaedlo cos. 

Cette enqu'ete de 1997 sur le revenu des m'enages a 'et'e pro duite par I’Office philippin da Statistique. 
courbe de Lorenz est pr'esent'ee en Figure 4.1. Fig. 4.1 - Courb e de Lorenzo.o 



Remarque4.1Surlegraphique,onindiquetoujoursladiagonale.LacourbedeLorenzest'egale'aladiagonalesitou 




##CourbedeLorenzetindicesd'irTegalit'e 

##Etapel:oninstallelapackageineq 

utils: ::menulnstallPkgs() 

#choisir'ineq'danslaliste 

##Etape2:onchargelepackageineq 

local({pkg<-select.list(sort(.packages(all.available=TRUE))) 

+ if(nchar(pkg))library(pkg,character.only=TRUE)}) 

#choisir'ineq'danslaliste 

##Utilisationdelabasededonn'eesllocos 

#Enqu y 'etesurlerevenuderOfficedeStatistiquePhilippin 

data(llocos)attach(llocos)#plot(Lc(income),xlab="proportiondemenages", 


ylab="proportionderevenu",main="")4.3.3lndicedeGiniL’indicedeGini,not'eGest'egaradeuxfoislasurf 


llestpossibledemontrerque:G=in(n-i)n/=in;=i|x/-xy|2'x.Enutilisantlastatistiqued'ordrex(i),...,x(/) / ... / x(n),l'indiced 


ecrire 

G=ln-12n/=i/x(/)n'x-(n+l).L’indicedeGiniestcomprisentre0etl.S'ilestprochede0,touslerevenussont'egaux.J 


l,lesrevenussonttr'esin'egaux.4.3.4lndicedeHooverL'indiced''equir'epartitiondeHoover(ouRo£>/n/-/c> 

'e'alacourbedeLorenz,carilestpossibledemontrerqu'ilcorrespond'alaplusgrandedistanceverticaleentrelacoui 


Lequintileshareratioestdefinitpar 


Q5R= |S2. 
520 

Ledecileshareratioestdefinitpar 


Cesquantit'essonttoujoursplusgrandesqueletaugmententavecl’iregalife' Cesdeuxrapportssont 
facilementinterpr'etables,parexemplteg3)SR=5, celasignifiequelerevenumoyendes20%lesplus 
richesest5foisplusgrandquelerevenumoyende20%lespluspauvres. 

4.3.6lndicedepauvret'eUnindicesimpledepauvret'econsiste'acalculerlepourcentagedelapopulation( 


moiti'edelam'ediane.4.3.7lndicesselonlespaySLetableau4.4reprendpourtouslespaysi'indicedeGini 


etlerapportdes20%lesplusrichessurles 

20%lespluspauvres.(r'ef'erence:UnitedNations2005DevelopmentProgrammeReport,page270). 


ExerciC6SExercice4.1'Etudiezlesproprret'es(circularit'e,r'eversibilit'e,identit' 


prixpr'esent'es.57 


eettransitivite)detouslesindicesde 



Tab.4.4-Mesuresdel'in'egalifeparpays 


7Slovakia25.8 6.7 

8BosniaandHerzegovina 26.2 5.4 

9Uzbekistan26.8 6.1 

10Finland26.9 5.6 

HHungary26.9 5.5 

12RepublicofMacedonia 28.2 6.8 

13Albania28.2 5.9 

14Germany28.3 6.9 

15Slovenia28.4 5.9 

16Rwanda28.9 5.8 

17Croatia29 7.3 

18Ukraine296.4 
19Austria307.6 
20Ethiopia306.6 
21Romania30.38.1 
22Mongolia30. 317. 89.1 
23Belarus30.46.94.6 
24Netherlands30.99.25.1 
25Russia317.14.8 
26SouthKorea31.67.84. 71998 
27Bangladesh31.86. 84.62000 
28Lithuania31.97.95. 12000 
29Bulgaria3 1.99. 95. 82001 
30Kazakhstan32. 37.55. 12003 
31Spain32.595.4199032lndia32.57.34.9199933Tajiki 


5.2200334France32.79.15.6199535Pakistan3; 


Chapitre5 

S'eriestemporelles,filtres,moyennes 
mobilesetd esaisonnalisation 


5.1Definitionsgeneralesetexemples 

5.1. ID 'efinitionsD'efinition5.1Unes'erietemporelleestunesuited'observationsd’unequantit' 

er epeteedansletemps. 

On'enonceeng'en'erairhypoth'esequelesintervallesdetempssont'equidBitifflfbfe5rl|H5relleestnffilE 

yi / ...,yt / ...,yT.Onnote'egalementT={l / 2 / ... / t,... / T>rensembledesinstantsauxquelslesobservationssontr' 

ealises. 

Lines 'erietemporellepeutsecomposerde:-unetendanceTt,-unecomposantecycliqueCt(nousn"etudieronspasc 
isoleretestimerunecomposantesaisonni'ere,etd'esaisonnaliserlas'erie,-r'ealiserunepr'evisionpourdesvaleu 
QTR:Quarter,trimestresdepuislelertrimestrel978jusqu'au4'emetrimestrel98559 



-DISH:Nombredelave-vaisselles(dishwashers)e>qides(milliers) 

-DISP:Nombredebroyeurd'ordures(disposers)expedes(milliers) 

-FRIG:Nombreder' efriGprateursexpdes(milliers) 

-WASH:Nombredemachines'alaver(washingmachine)e>piii , es(milliers) 

-DUR:D'epensesenbiensdurablesUSA(milliardsdedollarsdel982) 

-RES:lnvestissementresidentielprieUSA(milliardsdedollarsdel982) 


Tab.5.1-BiensmanufacturesauxUSA 


QTRDISH DISP 1 

FRIG 

WASH 

DUR 

RES 

1841 798 

1317 

1271 

252.6 

172.9 

2957837 

1615 

1295 

272.4 

179.8 

3999821 

1662 

1313 

270.9 

180.8 

4960858 

1295 

1150 

273.9 

178.6 

58948371271 


1289 

268.9 

174.6 

68518381555 


1245 

262.9 

172.4 

78638321639 


1270 

270.9 

170.6 

88788181238 


1103 

263.4 

165.7 

979286812771273 



260.6 

154.9 

1058962312581031 



231.9 

124.1 

1165766214171143 



242.7 

126.8 

1269982211851101248.6 



142.2 

1367587111961181258.7 



139.3 

1465279114101116248.4 



134.1 


1562875914171190255.5122.3 

165297349191125240.4110.4 

174807069431036247.7101.2 

1853058211751019249.1103.41955765912691047251.8100.120602837973918262.0115.8216588671102 

'erateursvendusamanifestementunecomposantesaisonni'ereetunetendance.EnlangageRQTR=c(l,2,3,4,5, 




1978 1980 1982 1984 1986 


Fig.5.1-D'epensesenbiensdurablesUSA(milliardsdedollarsdel982) 



1978 1980 1982 1984 1986 


Fig.5.2-Nombreder'efrig'erateursvendusdel978'al985 


DISP=c(798, 837, 821, 858, 837, 838, 832, 818, 868,623, 662, 822, 871, 791, 759, 734, 706, 

582,659,837,867,860,918,1017,1063,955,973,1096,1086,990,1028,1003) 

FRIG=C(1317, 1615, 1662, 1295, 1271, 1555, 1639, 1238, 1277, 1258, 1417, 1185, 1196, 1410, 1417, 919, 943, 1175 


1242, 1684, 1764, 1328)WASH=c(1271, 1295, 1313, 1150, 1289, 1245, 1270, 1103, 1273, 1031, 1143, 1101, 1181,1 
4, 161.9,159.9,170.5,173.1,170.3, 169. 6,170.3,172. 9,175,179. 4)plot(QTR,DUR,type="l")plot(QTR,FRIG,type= 





Tab.5.2-lndicedesprix'alaconsommation(France) 


pt 1970 

1971 

1972 

1973 

1974 

1975 

1976 

1977 

1978 

janvier 97.9 

102.5 

108.3 

115.5 

127.4 

145.9 

159.9 

174.3 

190.3 

f'evrier98.2 

103.0 

108.9 

115.8 

129.1 

147.0 

161.0 

175.5 

191.7 

mars98.5 

103.4 

109.4 

116.4 

130.6 

148.2 

162.4 

177.1 

193.4 

avril99.0 

104.0 

109.8 

117.2 

132.7 

149.5 

163.8 

179.4 

195.5 

mai99.4 

104.7 

110.4 

118.3 

134.3 

150.6 

164.9 

181.1 

197.4 

juin99. 8105.1 


111.0 

119.2 

135.8 

151.7 

165.6 

182.5 

198.9 

juilletl00.0105.6 


111.9 

120.2 

137.5 

152.8 

167.2 

184.1 

201.5 

acTutlOO. 4106.0 


112.5 

121.0 

138.6 

153.8 

168.4 

185.1 

202.5 

septembrel00.8106.5 

113.2 

122.1 

140.1 

155.1 

170.2 

186.7 

203.8 

octobrel01.2107. 1114.2 


123.4 

141.8 

156.3 

171.8 

188.2 

205.7 

novembrel01.6107.5114.9 

124.5 

143.1 

157.3 

173.2 

188.9 

206.8 

d 'ecembrel01.9 108.01 15. 5 125.3 

144.3 

158.2 

173.8 

189.4 

207.8 



Fig.5.5-Rapportenglissementannueldesindicesdeprixpi/pt-i2#lprix=c(97.9,98.2,98.5,99,99.4,99.8,100,100.‘ 






108.3,108.9,109.4,109.8,110.4,111,111.9,112.5,113.2,114.2,114.9,115.5, 

115.5.115.8.116.4.117.2.118.3.119.2.120.2.121.122.1.123.4.124.5.125.3, 

127.4.129.1.130.6.132.7.134.3.135.8.137.5.138.6.140.1.141.8.143.1.144.3, 
145.9,147,148.2,149.5,150.6,151.7,152.8,153.8,155.1,156.3,157.3,158.2, 
159.9,161,162.4,163.8,164.9,165.6,167.2,168.4,170.2,171.8,173.2,173.8, 

174.3.175.5.177.1.179.4.181.1.182.5.184.1.185.1.186.7.188.2.188.9.189.4, 
190.3,191.7,193.4,195.5,197.4,198.9,201.5,202.5,203.8,205.7,206.8,207.8) 
lprix<-ts(lprix,start=c(1970,l),frequency=12) 

plot(lprix)plot(lprix/lag(lprix,-l))plot(lprix/lag(lprix,-12))Exemple5.4Donn'eesdunombredevoyageurs-kilom 


etresendeuxemeclasseexpriraesenmillionsde 

kilom'etres.Tab.5.3-TraficdunombredevoyageursSNCF 


acfut sept. oct. nov. dec. 
2850 2090 1850 1630 2420 
2960 2190 1870 1770 2270 
2140 1870 1760 2360 
2160 1940 1750 2330 
2100 1920 1670 2520 
1936 1784 2391 
2085 1856 2553 


mois/ann'eejanv.rev.marsavrilmaijuinjuil. 

19631750156018202090191024103140 
19641710160018002120210024603200 
196516701640177021902020261031902860 
196618101640186019902110250030302900 
196718501590188022102110248028802670 
1968183417921860213821152485258126392038 
1969179818501981208521202491283427251932 
19701854182320052418221927222912277121532136 
19712008183521202304226421752928273821782137 
19722084203421522522231826842971275922672152 
197320812112227926612281292930892803229622102135 
197422232248242127102505302133273044260725252160 
197524812428259629232795328735983118287527542588 
1976266726682804280629763430370530532764280227073307 
1977270625862796297830533463364930952839296628633375 

1978282028573306333331413512374431792984295028963611197933132644287232673391368: 


1910 2537 
2009 2546 


1978 2723 
2862 


2876 

3266 



Fig.5.6-TraficdunombredevoyageursSIMCF63 




5.2Descriptiondelatendance 

5.2.1Lesprincipauxmodeles 

Plusieurstypesdemodlespeuvenfetreutiliaaspourdecrirelatendance. 

-Mod'elesd'ependantdutemps.LeEiedependdirectementdutemps.LemaaJepeutetreadditif: 
y t =f(t)+E t , 

oumultiplicatif y 

t=f(t)xE t . 

-Mod'elesexplicatifsstatiquesdasariechronologiqueelpenddesvaleursprisesparuneouplusieurs 

autress'erieschronologiques.y 

t =f(x t )+E t 

Lecaslin 'eaireestleplusfacile' atraiter yt=b 

o +bix t +E t- 

-Mod'elesauto-projectifs.Las'eriechronologiqueautempsfJdnddesespropresvaleurspaees 
yt= f(yt-i,yt-2,yt-3, ...,y t- P )+E t 

-Mod'elesexplicatifsdynamiques:las'eriechronologiqued'ependdesvalesHTB|b®setpaasesd'une 

oudeplusieursautress'erieschronologiques,parexemple:yt=)Li+0iyt-i+02yt-2+-"+0pyt-p+<j!)ixt-i+(p2Xt-2 

H \-(p qXf-q +E t. 

5.2.2Tendancelin 'eaireLatendancelaplussimpleestlin'eaire.Onpeutestimerlesparam'etresaurnoyenc 

ethodedes 

moindrescarr'es.C’estuner'egressionsimple.Tt=a+£>t.5.2.3TendancequadratiqueOnpeututiliserune 


desmoindrescarr'es.C'estuner'egressionavecdeuxvariablesexplicatives.Tt=a+bt+ct25.2.4Tendancepo 

'eaire,onpeutseramener'aunprobremelin'eaire.Enposantzt=l/Tt,onazt=l+ibe-atc64 



Zt+1 


l+be~ a(t+1) 


1 +be~ at e~ a 
c 

1+(1 +be~ at )e~ a -e~ a 


Enposanta= 

onobtientz 



cequiestunmod'eleauto-projectif.OnpeutalorsdlerminerlesvaleursdeaetjSparunesimpleiegression 

lin'eaire.Ensuiteond'eduitadelamani'eresuivante: 


a= -log/3, 

etcommea=l-e-ac=l-j8 

ond'eterminecparc=l-/3a.Enfin,onremarquequezt-lc=be-atc,onpeutd'eterminerautantdevaleursde£>quel' 



Fig.5.7-Exempledefonctionlogistiqueavecc=0.565 




5.30p' erateursdecfe'calageetdediffe'rence 

5.3.10perateursdecfecalage 

Afindesimplifierlanotation, onutilisedesoperateursdecfe'calageOnd efinitl'operateurdedecalage 
"retard"(enanglais/agoperator)Lpar 

Ly t =y t - 1 , 

etrop'erateur(enanglaisfon/i/ardoperator)"avance”F 


Fy t =yt+ 1 , 


I'op'erateuridentit'e/y 


t =yt. 


L'op'erateuravanceestl'inversedel'op'erateurretard 


FL=LF=I. 


Onpeutdonc 'ecrireF-i=LetL 

_1 =F. 

Ona'egalement-L2yt=/./.yt=yt-2,-Lc7yt=yt-c7,-F(jyt=yt+(?,-Lo=Fo=/,-/.-(?yt=F(7yt=yt+ £7 .5.3.20p"erateurdiff 
'eedanslegraphique5.9.EnlangageR##Tendancelin'eaireetdiff'erence#lin=10+0.3*(0:50) + rnorm(50,0,l)| 



0 10 20 30 40 50 


Fig .5 .8-S ' erieavecunetendancelieairedependantdutemps 
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Fig.5.9-Diff'erenced'ordreundelas'erieavecunetendanceleaire 


Onpeutconstruirerop'erateurdiff'erenced’ordredeuxerrelevantVaecarr' 

V2=VxV=/-2L+L2L'op'erateurdiff'erenced'ordredeuxpermetd'enleverunetendancequadratique.Eneffet,sila 

eries'ecrit 

yt=a+bxt+cxt2+Et,a\orsV2yt=(l-2L+L2)yt=a+bxt+cxt2+Et-2a-2bx(t-l)-2cx(t-l)2-2Et-i+a+bx(t-2) 

'erateurs,lacomposantesaisonni'eredispara''it.EnlangageR67 






1979 1980 1981 1982 1983 1984 1985 


Fig.5.10-Diff'erenced'ordre4delavariableventede'iefrigprateurs’ 


##Venteder'efrig'erateursdifrerenced'ordre4 

#FRIGm4=FRIG-lag(FRIG,-4)plot(FRIGm4)Exemple5.7Sionappliqueunediff'erencesaisonni'ered'ordrel2sur 


eesdunombredevoyageurs- 

kilom'etresytendeuxi'emeclasseexprim'eesenmillionsdekilom'etresdelaSfeltelndancesaisonieire 

dispara , 'it(voirFigure5.12).Onaainsilanouvellevariablezt=Vi2yt=(/-Li2)yt=yt-yt-i2.Uneautremani'eredefaire( 


erence, 

cequirevient' aprendrelelogarithmedurapportdelavariable(voirFigure5.13).Ond ' efinitainsi une 
nouvellevariablei/t:i/t=Vi 2 logyt=(/-Li 2 )logyt=logyt-logyt-i 2 =logytyt-i 2 .Time 



Fig. 5. H-TraficdunombredevoyageursSNCFEnlangageRtrafic=c(1750,1560, 1820, 2090, 1910, 2410,3140,28! 






1965 1970 1975 1980 

Time 

Fig.5.12-Diff'erenced'ordrel2surlaserietraficdunombredevoyageursSNCF 



1965 1970 1975 1980 

Fig.5.13-Logarithmedurapportd'ordrel2surlas'erietraficdunombredevoyageursSNCF 

1834,1792,1860,2138,2115,2485,2581,2639,2038,1936,1784,2391,1798,1850,1981,2085,2120,2491,2834 
806,2976,3430,3705,3053,2764,2802,2707,3307,2706,2586,2796,2978,3053,3463,3649,3095,2839,2966,: 
"eairesUnfiltrelin'eaired'ordrem=pi+p 2 estd'efiniparFL=p 2 y — P iWjL-j=w- P iL P i+w- P i+iL P i-i+---+w-iL+wol+' 





o'upi,p GNetw j ER. 


5.4.2Moyennesmobiles:definition 

Unemoyennemobiled'ordrem=p +p 2 + lestunfiltrelin'eairetelque 

pi 

wj = 1, pourtout/=,-p i 2. 

i — p i 

Beaucoupdemoyennesmobilesontdespoidj5ipfositifs,maispastoutes. 
Unemoyennemobileestsyetriquesipi =p 2 =p,e t 

1/14 =w-j, pourtouty=l,...,p. 
Unemoyennemobilesym'etriqueestditenon-paaBelsi 

n/y=cstpourtouty=,-p 2. 

5.4.3Moyennemobileetcomposantesaisonarrie 

Unemoyennemobileestunoutilint'eressantpourlisserunerfetemporelleetdoncpourenleverune 
composantesaisonni'ere.Onutilisedepr'ef'erencedesmoyennesmobilesnoei^HHiHordreegarala 
p'eriode,parexempled’ordre 7 pourdesdonn'eesjournali'eres,d’ordrel 2 pourdeedBrimBhsuelles.Par 
exemple,pourenleverlacomposantesaisonni'eredueaujourdelasemaine,onpeutappliquerunemoyenne 
mobilenon-pond'er'eed'ordre 7 .MM( 7 ) = 17 L 3 +L 2 +L+/+F+F 2 +F 3 


Cettemoyennemobileaccordelem , 'emepoids'achaquejourdelasemaine.Eneffet, 

MM( 7 )yt= 17 (yt- 3 +yt- 2 +yt-i+yt+yt+i+yt+ 2 +yt+ 3 ).Pourlescomposantessaisonni'eresd’unep'eriodepaire,ilr 


ees 

non-pond 'er'ees.llexistedeuxtypesdemoyennernobilecentr'ee:-Silap'eriodeestpaireet'egale' a m,(m= 4 poui 

mobiled'ordreimpairaccordantundemi-poidsauxdeuxextr'emit'es.Parexemple,pourdesdeBn' 

trimestrielles,lamoyennemobileestd'efinieparMM( 4 )= 18 L 2 + 2 L+ 2 /+ 2 F+F 2 .Ainsi,chaquetrimestreconserveler 

'eesseront''perdues''auxextr'emit'esdelas'eries.Exemple 5 . 8 LavariableVefrig'erateur"estliss'eegr''ace'a 



1978 1980 1982 1984 1986 

Time 

Fig.5.14-Nombreder' efrigprateursetmoyennemobiled'ordre4 


EnlangageRdec=decompose(FRIG)moving_average=dec$trendplot(FRIG)lines(moving_average)Unemoyer 


sonni'ere.5. 5Moyennesmobilesparticuli'eres5.5.iMoyennemobiledeVanHar 

'eriode4et5etconservelestendancespolynomialesjusqu''al'ordre3.5.5.3MoyennemobiledeHender 




MoyennemobiledeHendersond'ordre2m-3,o'um>4 

m + 1 

MM H = 6jU , 

j — m-1 

o' U0J = 

315((m-lf -i 2 )(m 2 ~/ 2 )((m+ 1) 2 -i 2 )(3rrf -16-11 i 2 ) 

8m(/7? -l)(4nf -l)(4nf -9)(4nf -25) 

MoyennemobiledeHendersond'ordre2m-3=5(m=4) 

1286(-2^L_ 84L+160/+84/ __ 21F 2) 

MoyennemobiledeHendersond'ordre2m-3=9(m=6) 

12431(-99L4-24L3-288L2+648L+805/+648F+288^ _ 24/r3 _ 99F 4) 

MoyennemobiledeHendersond'ordre2m-3 = ll(m=7) 

1923 78(— 2574L5— 2475L4+3300/-3+13050/-2+22050L+25676/ 
+22050F+13050F2+3300F3-2475F 4 -2574^) 

MoyennemobiledeHendersond'ordre2m-3 = 15(m=9)1193154(-2652L7-4732L6-2730L5+4641L4+16C 

+ 37422F+28182F2+16016F3+4641F4— 2730F5— 4732F6— 2652F 7 ) 

5.5.4M 'edianesmobilesSilesdonn'eescontiennentdesvaleursaberrantesouextr''emes,onpeutremplac 

unem'edianemobile.Parexemplelam'edianemobiled'ordre5estd'efiniepar:Med(5)t=M'ediane(yt-2,yt-i / yt,yt+ 

'eesetlatendanceSm=l/\-la(V'am-Tam).72 



Eng'eeral,onnedisposepasduramenombred'observations,pourchaquemois.On$rfecaunajustement 

afinquelasommedescomposantessaisoeressoitegale'aefo: 



Onpeutensuiteproceder'alacte'saisonnalisationdelaH'lepar 
Y am =Y am —5 m . 


5.6.2M 'ethodemultiplicative 

Soitunes'erietemporeller'egieparunneitelmultiplicatifdutype 

Yam =T am XS m XE am- 

o'ua=l,...,>4repr'esenteparexemplerann'eeet/T)=l,..,Mrepesenteparexemplelemois.Latendance 

estsuppos'eeconnuesoitparunajustement,soitparunemoyennemobile. 

Onisolelacomposantesaisonni'ereenfaisant,pourchaquemois,lamoyennedesrapportsentrelesvaleurs 

observ'eesetlatendance:Sm=L4-laY 

Tam' 

'Anouveau,onr'ealiseunajustementafinquelamoyennedescomposantessai§ft§v^oitegale'alOn 

corrigedonclescoefficientsSmparSm=SmliMmSm.Lad'esaisonnalisationser'ealisealorsparunedivisionV , am=V'amf 


erateursLeTa- 

bleau5.4contientlavariable'venteder'efrig'erateurs',lamoyennemobiled'ordre4,lacomposan4esaisonni' 
ets'eried'esaisonnalis'eeaumoyendelam'ethodeadditive.LeTableau5.6pr'esentelad'esaisonnalisationau 
moyendelam 'ethodemultiplicative. EnlangageRdeco=decompose(FRIG,type="multiplicative")plot(deco) 5 , 
'ediction'al'horizonl,etconsiste'ar'ealiserunemoyennedesvaleurspass'eesenaffectantdespoidsmoinsimpor 


Tab.5.4-D' ecompositiondelavariableFRIG,rrethodeadditiv< 


QTR 

FRIG 

MM 

FRIG-MM 

Desaison 

1 

1317 



1442.58 

2 

1615 



1505.13 

3 

1662 

1466.50 

195.50 

1451.20 

4 

1295 

1453.25 

-158.25 

1490.09 

5 

1271 

1442.88 

-171.88 

1396.58 

6 

1555 

1432.88 

122.13 

1445.13 

7 

1639 

1426.50 

212.50 

1428.20 

81238 

1390.13 

-152.13 

1433.09 


91277 1325.25 -48.25 1402.58 

101258 1290.88 -32.88 1148.13 

1114171274.13 142.88 1206.20 

1211851283.00 -98.00 1380.09 

1311961302.00 -106.00 1321.58 

1414101268.75 141.25 1300.13 

1514171203.88213.13 1206.20 

169191142.88- 223.88 1114.09 

179431095.00-152.00 1068.58 

1811751083.2591.75 1065.13 

1912691109.88159.13 1058.20 

209731150.88- 177.881168.09 
2111021218.50-116.501227.58 

2213441296.5047.501234.13 
2316411368.88272.131430.20 

2412251454.13-229.131420.092514291512.00-83.001554.582616991512.00187.001589.132717491475.] 





Tab.5.6-D' ecompositiondelavariableFRIG,methodemultiplicative 


QTR 

FRIG 

MM 

FRIG/MM 

Desaison 

1 

1317 



1453.85 

2 

1615 



1493.76 

3 

1662 

1466.50 

1.13 

1434.00 

4 

1295 

1453.25 

0.89 

1516.45 

5 

1271 

1442.88 

0.88 

1403.07 

6 

1555 

1432.88 

1.09 

1438.26 

7 

1639 

1426.50 

1.15 

1414.15 

81238 

1390.13 

0.89 

1449.70 

91277 

1325.25 

0.96 

1409.70 

101258 

1290.88 

0.97 

1163.56 

1114171274.13 

1.11 

1222.61 

1211851283.00 

0.92 

1387.64 


1311961302.00 0.92 1320.28 

1414101268.75 1.11 1304.15 

1514171203.881.18 1222.61 

169191142.880.80 1076.15 

179431095.000.86 1040.99 

1811751083.251.08 1086.79 

1912691109.881.14 1094.91 

209731150.880.851139.39 
2111021218.500.901216.51 
2213441296.501.041243.10 
2316411368.881.201415.88 

2412251454.130.841434.482514291512.000.951577.492616991512.001.121571.452717491475.131.191! 

'eeparunedroitehorizontale.Autrementdit,onsupposequeXr=a.75 




Fig.5.15-D' ecompositiondelaefiedeventesdereTrigprateurs5.1 

Decomposition of multiplicative time series 





Lelissageexponentielpeut'etreobtenuaumoyendelam'ethodedesmoindresefflinmnimisantena 

lecrit'ereQ=r-iy=o/3y(Xr-y-a)2.Enannulantlad'eriv'eeparrapport'aa,onobtient2r-iy=o/3y(XT-y-a)=0 / cequidonne 

'eeparunedroitequelconquedetypea+b(t-T).Onappliquealorsunlissageexponentieldoublepourobtenirlapr'e< 




Enannulantlesderiveespartiellesparrapport'aaetb,onobtient 


cequidonneQ 


1 ( Xr-j —a+bj) = 0 


T— 1 

2 0 (X T -j -a+bj)j=0. 


T ~ 1 . T “ 1 . . 

S ff X r-j ~a 0 +b jf3 J =0 


T-lj=OJ0 T- 1 T- 1 

X T -j -a jp j +b j 2 0 =0. 


i=oB j -- 


i= o/3 


»y=oy/3y=/3 

(1-/3) 2 

^oj2pj=p(i+m-PhonauuuununnunnT-y=opjXT-j-ai-p+bp(i-P)2=OT-ij=omxT-j-ap(i-i3)2+bp(i+m-^ 


EnnotantmaintenantSiTlas'erieliss'eeSiT=(l-/3)T-iy=o/3yXT-y,etS2Tlas'eriedoublementliss'eeS2T=(l-j3)T-iy=o/ 




Exemple5.1OLetableau5.8rendcompteduprixmoyendumazoutpourl00 (achatentre800etl500 
)enCHFpourchaquemoisde2004'a2007(Source:Officefecfe'raldelastatistique,2008). 


Tab.5.8-PrixmoyenduMazoutpourl00 (achatentre800etl500) 

mois/ann 'ee200420052006 2007 

f'evrier51. 5167. 3288. 7081.32 
mars55. 6075. 5288. 9282.06 
avril55. 7279.8392. 5888. 05 

mai58. 7173. 2293. 6588. 24juin58. 8275. 3891. 8888. 95juillet58. 4183. 9795. 3592. 10ao , 'ut64. 9284. 2395. 83? 
'ecurs\ve~Xt(l)=(l-l3)Xt+p''Xt-ia)\Xo(l)=Xi,o'uSit=''Xt(l).Onobt\ent:Sii='Xi(l)=(l-p)Xi+p~Xo(l)=(l 



Onobtient: 


Si =(l-j 8 )S l +/3Sq =(1— jB)S } +/3S! =S 1 1 =54.23, 

Sf =(1-/3 )S ^ +/SS? =0.3x53.414+0.7x54.23 = 53.99, 
S| =(1-JB)S 3 +/3S| =0.3x54.070+0.7x53.99 = 54.01, 


etainsidesuite.-Oncherchealors" 

X t (k)=a+bk 

pourchaquet.Onprendicik=l, X t (l)=a+bavec: 


a 

b 



°i- 3 s 1 -s 2 
0.7 ^ 


Letableau5.9rendcomptedesr'esultatspourlesaBBs2004'a2007. 

Lafigure5.16repr'esentelas'erieinitiale,lelissageexponentielsimpleetlelissageexponentieldoubleet 

peut"etreobtenueenlanguageRaumoyenducodesuivant: 


#Lissageexponentieldoubleaveck=lmazout=c(54.23,5 1.5 1,55.60,55.72, 58.7 1,58.82, + 58.41,64.92 ,63.95, 7 


+ 79. 83, 73. 22, 75. 38, 83. 97, 84. 23, 97. 29, 99. 31, 89. 88, 87. 18, +86. 16, 88. 70, 88. 92, 92. 58, 93. 65, 91. 88, 95. 35, 95. 
+p*liss2[i-l]#formuler'ecursive}#Lissageexponentieldoubleaveck=la=2*liss-liss2b=((l-p)/p)*(liss-liss2)lis 


Tab.5.9-Lissageexponentielsimpleetdoubledela£rietemporellePrixmoyenduMazoutpourl00litres 

(achatentre80 0etl500litres)enCHF 

~annee mois X t 5* Sf a b a+b 
2004 I 54.23 54.230 54.230 54.230 0.000 54.230 

2 51.51 53.414 53.985 52.843 -0.245 52.598 

3 55.60 54.070 54.011 54.129 0.025 54.154 

455.72 54.564 54.177 54.952 0.166 55.119 

558.71 55.808 54.666 56.950 0.489 57.440 

658.82 56.712 55.280 58.144 0.614 58.757 

758.41 57.221 55.862 58.580 0.582 59.163 

864.9259.531 56.963 62.099 1.101 63.199 

963.9560.857 58.131 63.582 1.168 64.750 

1072.9864.494 60.040 68.947 1.909 70.856 

1170.2566.22161.894 70.547 1.854 72.401 

1268.2466.82663.374 70.279 1.480 71.759 

2005163.0065.67864.065 67.292 0.691 67.983 

267.3266.17164.697 67.645 0.632 68.277 

375.5268.97665.98171.971 1.284 73.254 

479.8372.23267.85676.608 1.875 78.483 

573.2272.52869.25675.799 1.402 77.201 

675.3873.38470.49676.272 1.238 77.510 

783.9776.56072.31580.8051.819 82.624 

884.2378.86174.27983.4431.964 85.407 

997.2984.39077.31291.4673.033 94.501 

1099.3188.86780.77896.9533.466100.420 
1189.8889.17083.29695.0442.51897.562 
1287.1888.57384.87992.2671.58393.850 

2006186.1687.84985.77089.9280.89190.819288.7088.10486.47089.7380.70090.439388.! 

1.41380.486281.3284.03487.15680.911-1.33879.573382.0683.44186.04180.842-1.11479 
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temps 


Fig.5.16-EvolutionduprixdumazoutenCHF(achatentre800etl500),lissageexponentieldoubleetlissageexponer 


Exercices 

Exercice5.1DesaisonnalisezlaE3‘iesuivante(c'estuneerietrimestriellesur3ani®es) 

2417,1605,1221,1826,2367,1569,1176,1742,2804,1399,1063,1755 

parlam'ethodeadditive,enutilisantunemoyennemobilecTordre4. 

Exercice5.2EnlangageRutilisezlaserie"ldeaths"quiestunesriequisetrouvedanslepackagedebase 
"datasets". Usezladocumentation,puisdaisonnalisezcetteisieparlesmethodesadditiveetmultiplica- 
tive.82 



Chapitre6 


Calculdesprobabilfesetvariables 

al eatoirese . IProbabilit eS6.1.1 'Ev'enementUneexp'erienceestditeare 


esultatOnnotewunr' esultat 

possibledecetteexp'erienceareatoire.L’ensembledetouslesr'esultatspossible£fegkTefcemple,si 

onjettedeuxpi'ecesdemonnaie,onpeutobtenirlesr'esultatsQ={(P,P) / (FP) / (PF) / (FF)>,avecFpourface"etPpoui 


eriencea&atoire 

comme"avoirdeuxfoispile"ou"avoiraumoinsunefoispile".Formellement,unevenementestunsous- 

ensembledeQ.-L''ev'enement"avoirdeuxfoispile”estlesousensemble{(PP)>.-L"ev'enement"avoiraumoinsi 


L'ensembleQestappere'ev'enementcertain,etrensemblevide0estappere'ev'enementimpossible. 

6.1.20p'erationssurles 'ev'enementsSurles'ev'enements,onpeutappliquerlesop'erationshabit 
'enemenMest"obtenirunnombrepair"etr'ev'enementfi"obtenirunmultiplede3",i''ev'enement4nBestr'ev'< 


Exemple6.1L'experiencepeutconsister'ajeterunKjalors 

Q={1,2,3,4,5,6>, 

etun'ev'enement,adtest"obtenirunnombrepair".Onaaiors 

A={2A,6}etA={l,3,5}. 


6.1.3Relationsentrelesevenements 

'Ev'enementsmutuellementexclusifs 

Si>4n6=0onditque>4etSsontmutuellementexclusifs,cequisignifieque-4etSnepeuventpas 

seproduireensemble.Exemple6.2Sionjetteund'e,r'ev'enement"obtenirunnombrepair”etr' 

evenement"obtenirunnombre 

impair"nepeuventpas''etreobtenusenm , 'emetemps.llssontmutuellementexclCKd : Btrepart,si Ton 
jetteund'e,les'ev'enementS/4:"obtenirunnombrepair"n'estpasmutuellementexclusifaaAecilfement 
6:"obtenirunnombreinf'erieurou'egara3".Eneffet,rintersectionde>4et6estnon-videetconsisteen 
r'ev'enement"obtenir2''.lnclusionSi/4estinclusdans6,on'ecriMc:fi.Onditque/4implique6.Exemple6.3Sionje 


>4={2>ete={2,4 / 6;.Onditque^impiiquee.6.1.4Ensembledespartiesd'unensembleetsysten 

" eP(.)estuneapplicationdeAdans[0,l],telleque:-Pr(Q)=l,-Pourtoutensembled'enombrabled’'ev'enementsA 



Tab.6.1-Syst' emecompletcbvenements 



Apartirdesaxiomes,onpeutd'eduirelespropeitessuivantes: 

Propri'et'e6.1Pr(0)=O.D 'emonstrat7'onComme0estd'intersectionvideavec0,onaquePr(0u0)=Pr(0)+Pr(0 


Donc,Pr(0)=2Pr(0),cequiimpliquequePr(0)=O.* 


Propri'et'e6.2Pr(A)=l-Pr(A).D 'emonstrat/onOnsaitqueAuA=QetAnA=0.Ainsi,onaquePr(Q) = Pr(AuA)=Pr(y 


Propri'et'e6.3Pr(A)<Pr(S)s/Ace.D "emonstrat/onCommeAcfi,ona6=(6nA)uA.Maisonaque(SnA)nA=0.Ains 



D demonstration 
Onaque/AuB=Au(BnA), 

aveo4n(BnA=0). __ 

AinsiPr(>4uB)=Pr(/4)+Pr(Bn/4). 

Hreste'amontrerquePr(B)=Pr(Bn/A) + Pr(>AnB) 

MaisB=(BnA)u(Bn>A) 

avec(Bn/A)n(Bn/A) = 0 _ 

DoncPr(B)=Pr(BnA) + Pr(BnA) 

* 

Propri'et'e6.5Prn/=iyA/<n/=iPr(>A/)D 'emonstrationNotonsrespect\\/e(r\er\tBi=Ai,B2=(A2\Ai),B3=(A3\(AujA 

2)), 

B4 = (A4|(AlUA2U>A3)) Bn=(An\(AlUA2UA 3 U ■"UjA n- 1)). 

Comment L4/=n/=iB/,etqueB/nBypourtouty=/ / alorsPrn/=iB/=n/=iPr(B/). Deplus, comme,pourtout/,B/cA/,onaquePr(i 



6.1.6ProbabilitesconditionnellesetireEt3endance 

D'efinition6.3Soientdeu@\enementsAetB,siPr(B)>0,alors 


Pr(A\B)= 


Pr(/AnB) 
Pr (B) 


Exemple6.5Sionjetteund' e,etquel’onconsiderelesdeuxevenementssuivants: 
-/Alaprobabilit'ed'avoirunnombrepairet 
-Blaprobabilit'ed'avoirunnombresuperieurouegara4. 

Onadonc-Pr(>4)=Pr({2,4,6>)=12,-Pr(B) = Pr( {4 J 5,6>)=36=12,-Pr(y4nB) = Pr({4,6>)=26=13,-Pr(/4|B)=Pr(/\nB)P 


Pr(,4|B) = Pr(/4).OnpeutmontrerfacilementquesL4etBsontind'ependants,alorsPrt4nB)=Pr(/4)Pr(B).6.1.7Th e 


Pr(B)=n;=iPr(/4;)Pr(B|/4;).Tab.6.2-lllustrationduth'eor'emedesprobabilit'estotales 



Eneffet,n/=iPr(A)Pr(B|A)=n;=iPr(BnA).Commeles'ev'enementsAnBsontmutuellementexclusifs,n;=iPr(Bn/4;)= 




Th 'eaxme6.2(deBayes)SoitA n unsysl£mecomplet(£\enements,alors 


Pr(>4/|B) = 


Pr(Ai)Pr(B\Ai) 

U PnAj)Pr(B\Aj)- 


Eneffet,parleth'eoe'medesprobabil&stotales, 


Pr(Aj)Pr(B\Aj) 


Pr(BnA ,) 
Pr(B) 


= Pr(A / |B). 


Exemple6.6Supposonsqu'unepopulationd'adultessoitcompeede30%defumeurs(>Ai)etde70% 
denon-fumeur(/A 2 ).NotonsBr'ev'enement"mourird'uncancerdupoumon"Supposonsenoutreque 
laprobabilit'edemourird'uncancerdupoumonest'egale'aPr(B|,4i) =20%si I’onestfumeuretde 
Pr(B|>A 2 )=l%sironestnon-fumeur.Leth' eoemedeBayespermetdecalculerlesprobabil&sapriori, 
c'est-'a-direlaprobabilit'ed'avoir'et'efumeursionestmortd’uncancerdupoumon.Eneffet,cette|Erobabilit' 
estnot'eePr(/4i|B)etpeut , 'etrecalcul 'eeparPr(/4i|B)=Pr(,4i)Pr(B|/4i)Pr(/4i)Pr(B|v4i)+Pr(,42)Pr(B|/42)= 

0^ = 0-06 

0.3x0.2 + 0.7x0.01 0.06+0.007 

Laprobabilit'edenepasavoir'et'enonfumeursionestmortd'uncancerdupoumonvautquant'aelle: 

Pr(A2\B)=P r{A2)Pr(B\A2)Pr(Ai)Pr(B\Ai)+Pr(A 2)P r{B\A2)=0.7x0.01 = 0.07 

0.&l).2+0.|tl0.01 0.06+0.007 0 ’ 


6.2Analysecombinatoire6.2.iintroductionLa nalysecombinatoireestl' 'etudemath 'emat 


natoireestunoutilutilis'edanslecalculdesprobabilit'es.6.2.2Permutations(sansr'ep'etition)Une|: 


rensemble{l / 2,3>.llexiste6mani'eresd'ordonnercestroischiffres:{l,2,3> / {l / 3,2>,{2,l,3>,{2 / 3 / l>,{3 / l / 2>,{2 



6. 2. 3Permutationsavecie petition 

Onpeutegalementseposerlaquestiondunombredemenasderangerdesobjetsquinesontpastous 

distincts.Supposonsquenousayons2boulesrouges(eetR)et3boulesblanches(nofeesS).llexistelO 

permutationspossiblesquisont: 

{R,R,B,B,B}, {R,B,R,B,B}, {R,B,B,R,B},{R,B,B,B,R}, {B,R,R,B,B}, 

{B,R,B,R,B},{B,R,B,B,R}, {B,B,R,R,B},{B,B,R,B,R}, {B,B,B,R,R}. 

Sil’ondisposedenobjetsappartenant'adeuxgroupesdetailleanetn 2 ,lenombredepermutationsavec 
r'ep'etitionestn! 


mW 

Parexemplesirona3boulesblancheset2boulesrouges,onobtient 

n!m!n 2!=5! _ 120 _ 

2 ! 3 ! 2X6 

Sirondisposedenobjetsappartenant'apgroupesdetailles/an 2 ,...,n p ,lenombredepermutations 
avecr'ep'etitionestn!ni!n 2 !x---xnp!. 


6.2.4Arrangements(sansr"ep"etition)Soitnobjetsdistincts.Onappelleunarrangementunemani'e 

electionnertobjetsparmilesnet 

delesrangerdansdesbo^itesnum'erot'eesdel'a/c.Danslapremi'erebcMte.onpeutmettrechacundesnobjets.Dc 

ite, onpeutmettre 

chacundesn-lobjetsrestants,danslatroisi'emebo , 'ite,onpeutmettrechacundesn-2objetsrestantset 

ainsidesuite.Lenombred’arrangementspossiblesestdonc'egara:/Afcfi=nx(n-l)x(n-2)x-"X(n-/f+l)=n!(n-/c-; 


tenircomptedeleurordre.Lenombredecombinaisonsestlenombredesous-ensemblesdetaille/cdansun 

ensembledetaillen.Soitrensemble{l,2 / 3,4,5>.llexistel0sous-ensemblesdetaille3quisont:{l / 2,3>,{l,2 / 4>,{'l, 

'ements,sil'onneveutpartenircomptederordreondivise>Afcnpar/f!.Lenombredecombinaisonsdetobjetsparminvi 


6.3Variablesal eatoires 

6.3.1Definition 

Lanotiondevariableafe'atoireformalisel'associationd’unevaleurffiaciltatd'uneexpfienceaiatoire. 

D'efinitior\6.5UnevariableslatoireXestuneapplicationdel'ensemblefondamen£tdlansR. 

Exemple6.70nconsid'ereuneeKpiencealatoireconsistant'alancerdeugpesdemonnaie.L'ensemble 

desr'esultatspossiblesestQ={(FF),(FP),(P,F),(P,P)>. 

Chacundes'erementsdeQauneprobateili^.Unevariablealeatoirevaassocierunevaleur'achacundes 

'erementsdeQ.Consid'eronslavariabtealirerepesentantlenombrede"Faces"obtenus: 

X=0 

□□0 avecuneprobabilttl/4 
lavecuneprobabilit'el/2 
2avecuneprobabilit'el/4. 

C'estunevariableaTeatoirediscr' etedontladistributiondeprobfflb&ittpresen$eenFigure6.:i. 


o 


Fig.6i-Distributionde"faces"obtenus.6.4Variablesareatoiresdiscr'etes6.4.iD'efini 

(x-ju)2=xezpx(x)x2-^2.Onpeutaussicalculerlesmomentsettouslesautresparam'etres.90 


6.4.2Variableindicatriceoubernoullienne 


LavariableindicatriceXdeparametrep6[0,l]aladistributiondeprobabilitessuivante: 

lavecuneprobabilifep 

Oavecuneprobabilifel-p. 

L'esp'erancevautp=E(X)=Ox(l-p) + lxp=p, 
etlavariancevauta 2 =var(X)=E(X-p) 

2 =(l-p)(0-p) 2 +p(l-p) 2 =p(l-p). 

Exemple6.80ntireauhasardunebouledansuneurnecontenantl8boulesrougesetl2boulesblanches.Si 

Xvautlsilabouleestrougeet0sinon,alorsXauneloibernoulliennedeparametrep=18/(18+12)=0.6. 


6.4.3VariablebinomialeLavariableareatoirebinomialedeparam'etresnetpcorrespond'arexp' 

eriencesuivant£)nrenouvelle 

nfoisdemani'ereind'ependanteune'epreuvedeBernoullidepatiaap,'o'upestlaprobabilitedesucces 
pouruneexp'erience'erementaire.Ensuite,onnoteXlenombredesob<tenus.Lenombredesucesest 
unevariableal 'eatoireprenantdesvaleursenti' eresdeO' anetayantunedistributionbinomiale. 
UnevariableXsuituneloibinomialedeparam'etreO<p<letd'exposantn,si 

?r(X =x) = nxpxqn-x,x= 0, l,...,n-l,n, 

o'ug=l-p,etnx=n!x!(n-x)!.Demani'eresynth'etique,siXaunedistributionbinomiale,onnote:>>C~8(n,p).Rapp< 


dubin / 'omedeNewton(p+g)n.(p+q)o=l(p+g)i=p+g=l(p+g)2=p2+2pg+g2=l(p+g)3=p3+3p2g+3pg2+g3=l(/ 


L’esperancesecalculedelamararesuivante: 


E(X) = xPr(X=x) 

x=0 

= n x " 


x /fcf x (onpeutenleverletermex=0) 

: i x 

n- 1 


x-1 

=npnx=m-l 


If cl' 


-1 rv 

=np n -iz=ori-lzp 


lJn-l)-(x-l) 


z cf n 1( z (enposantz=x-l) 
=np(p+g)n-i=np.Lavarianceestdonn'ee(sansd'emonstration)parvar(X)=npg.Exemple6.90ntireauhasarda' 


18boulesrougesetl2boulesblanches.SiXestlenombredeboulesrougesobtenues,alorsXauneloi 
binomialedeparam'etrep=18/(18+12) = 0.6,etd'exposantn=5.Donc,Pr(X=x)=5x0.6x0.45-x,x=0,l,...,4,5,cequ 
'eedanslaFigure6.2.Exemple6.10Supposonsque,dansunepopulationd''electeurs,60%des'electeurss'appr''( 



« 1 
° -\ 


0 1 2 3 4 5 

Fig.6.2-Distributiond'unevariableal eatoirebinomialeavecn=5etp=0.6. 

6.4.4VariabledePoissonLavariableXsuituneloidePoisson,deparam' 

etre^GR + si 

p r ( x=x ) =e-^ i=0 ,l,2,3 

OnnotealorsX~P(l).Lasommedesprobabilit'esestbien egale'al,eneffet 

oo X =oPr(X=x)=oox=oe-/\/\xx! =e-A» 0 1 x 

— =e“V = 1. 

x! 

L'esp'eranceetlavarianced'uneloidePoissonsont'egalesauparam'etre^.Eneffet 

E(X)=oox=oxPr(X=x)=<x»x=oxe-^xx!=e-/\<»x=ix^xx!=e-^oox=i/Ax-i(x-l)!=e-/\^ooz=oAzz!enposantz=x-: 


=e-AAe 4 



Fig.6.3-Distributiond'unevariabledePoissonavec^ = l. 


lwd=3,xlab="",ylab="",main="",frame=FALSE) 
#PoissonP(l)plot(dpois(0:7,l),type="h",lwd=3,xlab="",ylab="",main = "",frame= FALSE) 


6. 5Variableareatoirecontinue6. 5. ID definition, esp 'eranceetvarianceunevar 


Laprobabilit'equ'unevariableareatoirecontinuesoitinf'erieure'aunevaleurpBgniesditibnneepar 

safonctionder'epartition.PrP(<x)=F(x).Lafonctionder'epartitiond'unevariableareatoirecontinueesttoujours: 

'erieure'aunevaleurquelconquevaut:Pr[X<a]=a-oo/ : (x)c/x=F(a).DanslaFigure6.4,laprobabilit'ePrtX<a]estl'aire 



Fig.6.4-Probabilit' equelavariablealeatoiresoitinfe'rieure'aa 

Silavariableareatoireestcontinue,laprobaleid|iiti'elleprenneexactementunevaleurquelconqueestnulle: 

Pr[X=a]=0. 

L’esp'eranced'unevariableareatoirecontinuelaBtdpar: 

E(X)=. 

-oo xf(x)dx, 


etlavariancevar(X)=oo-oo(x-ju) 2 f(x)dx. 


6.5.2VariableuniformeUnevariableareatoireXestditeuniformedansunintervalle[a,b],(aveca<b)sisar' 

epartitionest: 

/ r (x)=DDDOsix<a(x-a)/(b-a)sia<x<blsix>b.Sadensit'eestalorsf(x)=DDDOsix<al/(ib-a)sia<x<ibOsix>b.Onpe 


Resultat6.2 


a 2 =var(X)= 


(b-a) 2 
12 ' 


D 'emonstrationDemani' ereg 'en 'erale,unevariancepeuttoujourss' ' 

ecrirecommeunmoment'arorigined'ordre2moinsle 


carr'edelamoyenne.Eneffet, 

a 2 =var(X) 

~ b a(x-,u) 2 f(x)dx 


=ba(X2+fJ 


2 -2x/j)f(x)dx 


=baX2 f(X ) dx + 


H 2 f(x)dx-2n xf(x)dx 


=baX2f(x)dX+JL/2-2jU 


=baX2f(x)dx-jL/2.0ncalculeensuiteunmoment'arorigined’ordre2:baX2f(x)cfx=baX2lb-ac/x=l£)-afaaX2dx=ljb-a; 



f(x) 



a b 

Fig.6.5-Fonctiondedensit'ed'unevariableuniforme 

a.b & 

Fix) 

* 

Fig.6.6-Fonctionder'epartitiond'unevariableuniforme 

6.5.3VariablenormaleUnevariableareatoireA'estditenormalesisadensit'evautf f j,c72(x)=laV27Texp-12x 


o'ujL/ERetaGR 2 sontlesparam'etresdeladistribution.Leparam'etre/jestapperelamoyenneetle 

param'etreal"ecart-typedelavariablenormale.p-oo+o 0 ^+ap-aFig.6.7-Fonctiondedensit'ed’unevariablenorm< 




Lafonctionderepartitionvaut 


F I J,C7 2 M — 


o 


du. 



-00 + 00 p-CT P jU + CT * 

Fig.6.8-Fonctionder' epartitiond'unevariablenormale 

6.5.4Variablenormalecentr"eer "eduiteLavariableareatoirenormalecentr'eer'eduiteestunevari 

d’esprancenullejL/=Oetde 

varianceCT2=l.Safonctiondedensit'evaut/b,i(x)=lV2nexp-X22.etsar'epartitionvaut<D(x)=Fo,i(x)=x-<x»lV2nexp 



Enposant 



onobtientu=za+ju,etdoncdu=adz.Donc, 


F ij,o 2 {x)— 


z 1 2 

2 


adz = O 


a 


Lestablesdelavariablenormalenesontdonneesquepourlanormalecentreereduite.Lestablesne 
donnent<D(x)quepourlesvaleurspositivesdex, carlesvaleursnegativespeuvenifretrouveesparla 
relationdesym etrie. 6.5.5 Distri but ionexponentielle 


UnevariableareatoireXaunedistributionexponentiellesisafonctiondedereafetdonneepar: 

f(x)=^exp-(^x), six >0 
Osinon 

Leparam'etre^estpositif.Quandx>0,safonctionder'epartitionvaut: 

F(x)=xof(u)du=xo\e-\udu=-e-xu x . 

= i _p - / ' x 
o j. t? 

Onpeutalorscalculerlamoyenne:R'esultat6.4E(X)=l/^D"emonstrationE(>X')=ooox/ : (x)dx=ooox^e-Axdx=-l+x2 


1 

A' 


* 

llest'egalementpossibledemontrerquelavariancevaut:var(X)=l/\ 2 . 6 . 6 Distri blltionbivari eeDe 



Fig.6.10-Fonctiondedensit'ed'unevariableexponentielleavec>A=l. 

Aveclesdistributionsmarginales,onpeutd'efinirlesmoyennesmarginales,etlesvariancesmarginales: 

/jx= oo- ooxfx(x) dx, et/jy= ■»- <»y fy(y) dy, 

02x=™-°o(x-iux)2fx(x)dx,et<j2Y=oo-oo(y-ijY)2fY(y)dy. 

Onappelledensit'esconditionnellesJesfonctionsf(x|y)=f(x,y)/Y(y)etf(y|x)=f(x,y)/x(x).Aveclesdistributionscond 


tionnelles:)Lix(y)=oo-<x,xf(x|y)c/x,etjuy'(x)=<x»-ooyf(y|x)c/y,CT 2 x(y)=oo-oo-fx-/jx(y)> 2 f(x|y)dx,eta 2 y'(x)=oo- 0 o-fy- J uy(x)> 2 f( 

'ejointedesdeuxvariables,alorsXetY'sontind'ependantssi/xY(x,y)=/ : x(x)fr(y) / x,y€=R-100 




6.7Propri etesdesesperancesetdesvariances 

DemanieregeneralepourdesvariablesafeatoiresXetr, etavecaetbconstants, onalesr'esultats 
suivants.R'esultat6.5E(a+bX)=a+bE(X) 


D'emonstrationE(a+bX)=R(a+bx)f(x)dx=a 


f(x)dx+b xf(x)dx=a+bE(X). 

R R 

R'esu\tat6.6E(aY+bX)=aE(Y)+bE(X) 

D'emonstrationE(aY+bX)=RR(ay+bx)f(x,y)dxdy 


=aRRyf(x,y)dxdy+bRR xf(x,y)dxdy 

=aRyRf(x,y)dxdy+bRXRf(x,y)dydx 

=aRyf(y)dy+bRxf(x)dx=aE(Y)+bE(X)* 


R'esultat6.7var(a+bX)=b2var(X).D'emonstrationvar(a+bX)=R[a+bx-E(a+bX)]2f(x)dx=R[a+jbx-(a+bE(> 



Demonstration 


var(X+r) = [x+y-E(X+Y)] 2 f(x,y)dxdy 

R R 

= [x-E(X)+y-E(Y)] 2 f(x,y)dxdy 

R R 

= R [x-E(X)] 2 +[y-E(r)] 2 +2[x-E(X)][y-E(Y)]f(x,y)dxdy 

R 

= va r(X) + va r( Y) + 2 cov(X, Y) 


Resultat6.9 Deplus, siXetYsontind'ependantes,onaf(x,y)=f(x)f(y)pourtoutx,y 
E(XY)=E(X)E(Y). 

D ' e m o n st rat i o n E (XY) = RRxyf(x) f(y)dxdy 


* 


= Rxf(x)dxRyf(y)dy 
=E(X)E(Y).* 

Enfin,siXetY'sontind'ependantes,onacov(X,y)=0,etdoncvar(X+Y')=var(X)+var(Y').Enfin,ilestpossibledecalc 


ependantes, 

et\dent\que(r\er\td\stribu'ees.Th'eor'eme6.3SoitXi,...,Xnunesuitedevariablesareatoires,ind'ependanteseti 

ees 

etdontlamoyennepetlavariancea 2 existeritetsontfinies,alorssi~X=lnni=iXi,onaE(~X)=p,etyar(~X)= 02 n.D'emo 


6.8Autresvariablesaleatoires 

6.8.1Variablekhi-carnee 

Soitunesuitedevariablesadatoiresindpendantes,normales,ceaT4Efe'duites,Xi,...,X p ,(c'est-'a-dire 
demoyennenulleetdevariancegale'al),alorslavariableadatoire 


p 



estappereevariableareatoirekhiecaprdegesdeliberte. 

llestpossibledemontrerqueE(x 

2 P )=P. 

etquevar(x 

p)=2 P- 



Fig.6.11-Densit'ed'unevariabledechi-carr'eavecp=l,2,...,106.8.2VariabledeStudentSoitunevariable 


es 

delibert'e,ind'ependantedeX,alorslavariableareatoiretp=Xx 2 p/pestappereevariableareatoiredeStudent'a/: 



Fig.6.12-Densit'esdevariablesdeStudentavecp=l,2et3etd'unevariablenornnalel03 



6.8.3VariabledeFisher 


Soientdeuxvariableseb'toireskhi-caesmdependante^'q-fespectivement'apetqide^deliberlB, 

alorslavariablealeatoire 

F - XpP 

Fp ' q ~xVP 


estappereevariablffiBdtoiredeFisher'apetqdegiesdeliberte. 



Fig.6.13-Densit'ed'unevariabledeFisher 

Remarque6.1llestfaciledemontrerquelecarr'ed’unevariabledeStuden1:'agd93jideliberteestune 
variabledeFisher'aletgdegr'esdelibert'e.6.8.4Variablenormalemultivari "eeLevecteurdevariable 


eedemoyenne 

/Li=(pi / ...,pp)etdematricevariance-covarianceZ(onsupposeparsimplicit'equeZestdepleinrang), 
sisafonctiondedensit'eestdonn'eeparfx(x) = l(2n)p/2|Z|i/2exp-12(x-^i)Z-i(x-/Li) / (6.2) 


pourtoutxeRp.Fig.6.14-Densit'ed'unenormalebivari'eeRemarque6.2Sip=l,onretrouvel'expression(6.1).10 


Uncas particulierest important: supposons)uelamatricevariance-covariancepeut'ecrireZ= 
diag(a2i,...,c|) / cequisignifiequetouteslescomposantesduvecteurXsontnon-coarlees.Danscecas, 


(2tt)p* \ 2\ 1/2 



(2 U of)™ 


exp- -(x-ju)Z (x-jli) 


□ □ 

1 - expD- P {XJ ~£ )2 a 


(2 nr*( U°>) j=i 


(2n)p/2( f =1 Oj ) j=1 exp 2g* 
= P y=il(27T)i/2ayexp (x . _ jJ . p 

^T“ 

=pj = 1 fXj(Xj), 0 ' u fxj(xj) = 1(2 TT02j) l/2exp - (Xj—fJ 


(Xj -m) 2 




2a2 

estladensit'edelavariableXy.Onconstateques'ilyaabsencedecone'lationentrelesvariablesnormales, 

alorsladensit'eduvecteurnormalpeuts’'ecrirecommeunproduitdede9 , i5[fianslecasmultinormal(et 

seulementdanscecas),rabsencedecorr'elationimpliquecloncrind'ependancedesvaeiatelBrear 

Demani' ereg 'en 'erale,siXestunvecteurdevariablesal 'eatoiresdemoyennejuetdematricevariance- 
covarianceZ,etsiAestunematricegxpdeconstantes,alorsE(AX)=AE(X)=Aji,etvar(AX)=Avar(X)A=AZA.Dansl 


pendantsamatricevariance-covariancen'estpasn 'ecessairementdepleinrang). Donc,siXestunvecteurmultinor 



Exercice6.2DeterminezlesvaleursydelavariablenormalecenfeiereduiteZtellesque: 


1. Pr[Z</]= 0,9332; 

2. Pr[-y<Z<y']=0,3438; 

3. Pr[Z<y]= 0,0125; 

4. Pr[Z>y]= 0,0125; 

5. Pr[/<Z<3]=0,7907. 


Exercice6.3SoitunevariablealeatoireX~A/(53;a 2 =100)representantleiesultatd’unexamenpour 
un'etudiantd'unesection.D'eterminezlaprobadpbbrqueler'esultatsoitcomprisentre33,4et72,6. 


Exercice6.4Soitunevariableal 'eatoire7C/V(50;o 2 =100). Determinezlepremierquartiledecettedis- 
tribution.Exercice6.5Ensupposantquelestaillesencmdes'etudiantsd'unpaysadmettentladistributionnormale 


A/(172;a2=9).Ondemandeded 'eterminerlepourcentageth 'eorique: 

a)d"etudiantsmesurantaumoinsl80cm.b)d"etudiantsdontlatailleestcompriseentrel68etl80.Exercice6.6S 


elavitesse 

detouteslesautomobilespendantunejoum'ee.Ensupposantquelesvitessesrecueilliessoientdisie&ir 

normalementavecunemoyennede72km/hetun'ecart-typede8km/h,quelleestapproximativementla 

proportiond'automobilesayantcommisunexc'esdevitesse?Exercice6.7Pourrassemblaged’unemachine,onpt 


esune 

loinormaledemoyennelOcmetd''ecart-typeO,2cm.Ongroupelescylindresen3cat'egories:A:d'efectueuxetinut 

'ee'al7degr'esdelibert'e;3.d'unevariabledeStudent'a8degr'esdelibert'e;4.d’unevariabledeFisher(uniquerr 


Chapitre7 

Tablesstatistiques 


Tab.7.1-Tabledesquantilesd'unevariablenormalecenteereduite 



Tab.7.2-Fonctionder' epartitiondelaloinormalecenteereduite 
(Probabil&detrouverunevaleurinlerieur'au) 



u 0 . 
0 . 0 . 
0 . 1 . 
0 . 2 . 
0.3. 
0.4. 


;0.01.02.03 .04 

:5000. 5040. 5080. 5120. 5160 

.5398.5438.5478.5517.5557 

.5793.5832.5871.5910.5948 

.6179.6217.6255.6293.6331.6368 

.6554.6591.6628.6664.6700.6736 


.05 .06 .07 .08 .09 


.5199 

.5596 

.5987 


.5239 

.5636 

.6026 

.6406 

.6772 


.5279 

.5675 

.6064 

.6443 

.6808 


.5319 

.5714 

.6103 

.6480 

.6844 


.5359 

.5753 

.6141 

.6517 

.6879 


0.5.6915.6950.6985.7019.7054.7088.7123 

0.6.7257.7291.7324.7357.7389.7422.7454 

0.7.7580.7611.7642.7673.7704.7734.7764 

0.8.7881.7910.7939.7967.7995.8023.8051.8078 

0.9.8159.8186.8212.8238.8264.8289.8315.8340 


.7157 

.7190 

.7224 

.7486 

.7517 

.7549 

.7794 

.7823 

.7852 


.8106 

.8133 


.8365 

.8389 


1.0.8413.8438.8461.8485.8508.8531.8554.8577 .8599 .8621 

1.1.8643.8665.8686.8708.8729.8749.8770.8790.8810 .8830 

1.2.8849.8869.8888.8907.8925.8944.8962.8980.8997 .9015 

1.3.9032.9049.9066.9082.9099.9115.9131.9147.9162 .9177 

1.4.9192.9207.9222.9236.9251.9265.9279.9292.9306.9319 


1.5.9332.9345.9357.9370.9382.9394.9406.9418.9429.9441 

1.6.9452.9463.9474.9484.9495.9505.9515.9525.9535.95451.7.9554.9564.9573.9582.9591.959! 

43.9945.9946.9948.9949.9951.99522.6.9953.9955.9956.9957.9959.9960.9961.9962.9963.9964 

7.9997.9998108 
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Tab.7.4-Tabledesquantilesd'unevariablex 2 'andegesdeliberte 


0.01 


ordreduquantile 
0.025 0.05 0.95 


0.975 0.99 


n = 1 I 
2 


0.000157 0.000982 0.003932 3.841 
0.02010 0.05064 0.103 5.991 

0.216 
0.484 
0.831 
1.237 
1.690 
2.180 


0.115 

1297 

1.554 

1.872 

1.239 

1.646 

2.0882.700 


io: 
11 : . 
122 . 
13 ^, 
14 ^, 
152 . 
165 . 
lit. 
18 


192 


.5583.247 

.0533.816 

.5714.4045.226 

.1075.0095.892 

.6605.6296.571 

.2296.2627.261 

.8126.9087.96226.30 

.4087.5648.67227.59 

.0158.2319.39028.87 

.6338.90710.1230.14 


5.024 6.635 

7.378 9.210 


0.352 7.815 

0.711 9.488 


9.348 11.34 

11.14 13.28 


1.145 11.07 

1.635 12.59 


2.167 14.07 

2.733 15.51 


12.83 15.09 

14.45 16.81 


16.01 18.48 

17.53 20.09 


3.325 16.92 19.02 21.67 


3.940 

4.575 


18.31 20.48 23.21 

19.68 21.92 24.72 

21.03 23.34 26.22 

22.36 24.74 27.69 

23.68 26.12 29.14 

25.00 27.49 30.58 

28.85 32.00 

30.19 33.41 

31.53 34.81 

32.85 36.19 


20t . 
211 ;. 
22S . 

23] l 

24] I 
251.. 
1442! 


37.57 

38.93 


.2609.59110.8531.4134.17 
.89710.2811.5932.6735.48 
.54210.9812.3433.9236.7840.29 
0.2011.6913.0935.1738.0841.64 
.0.8612.4013.8536.4239.3642.98 

.5213.1214.6137.6540.6544.312612.2013.8415.3838.8941.9245.642712.8814.5716.1540 
!5. 1527. 5729. 7960. 4864. 2068. 714626. 6629. 1631. 4462. 8366. 6271. 204828. 1830. 7533. 106! 



Tab.7.5-Tabledesquantilesd'unevariabledeStudent'andegesdeliberte 

ordreduquantile 
0.95 0.975 0.99 0.995 

n = 1 6.314 12.71 31.82 63.66 

2 2.920 4.303 6.965 9.925 

3 2.353 3.182 4.541 5.841 

4 2.132 2.776 3.747 4.604 

5 2.015 2.571 3.365 4.032 

6 1.943 2.447 3.143 3.707 

7 1.895 2.365 2.998 3.499 

81.860 2.306 2.896 3.355 

91.833 2.262 2.821 3.250 

101. 8122. 228 2.764 3.169 

111.7962.201 2.718 3.106 

121. 7822. 179 2.681 3.055 

131. 7712. 1602. 650 3.012 

141.7612.1452.624 2.977 

151. 7532. 1312. 602 2.947 

161. 7462. 1202. 5832.921 
171. 7402. 1102. 5672. 898 
181. 7342. 1012. 5522. 878 

191.7292.0932.5392.861201.7252.0862.5282.845211.7212.0802.5182.831221.7] 

.3812.648801.6641.9902.3742.639901.6621.9872.3682.6321001.6601.9842.3642 


Tab.7.6-Tabledesquantilesd’ordre0.95d'unevariabledeFisher'anL etn 2 degesdeliberte 


m= 12345 6 7 8 9 10 12 14 16 20 30 oo 

02 = 1161.4199. 5215. 7224. 6230. 2234. 0236. 8238. 9240. 5241. 9243. 9245. 4246. 5248. 0250. 1254. 3 
218.5119.0019.1619.2519.3019.3319.3519.3719.3819.4019.4119.4219.4319.4519.4619.50 
310.139.5529.2779.1179.0138.9418.8878.8458.8128.7868.7458.7158.6928.6608.6178.526 
47.7096.9446.5916.3886.2566.1636.0946.0415.9995.9645.9125.8735.8445.8035.7465.628 
56.6085.7865.4095.1925.0504.9504.8764.8184.7724.7354.6784.6364.6044.5584.4964.365 
65.9875.1434.7574.5344.3874.2844.2074.1474.0994.0604.0003.9563.9223.8743.8083.669 
75.5914.7374.3474.1203.9723.8663.7873.7263.6773.6373.5753.5293.4943.4453.3763.230 
85.3184.4594.0663.8383.6873.5813.5003.4383.3883.3473.2843.2373.2023.1503.0792.928 
95.1174.2563.8633.6333.4823.3743.2933.2303.1793.1373.0733.0252.9892.9362.8642.707 

104.9654.1033.7083.4783.3263.2173.1353.0723.0202.9782.9132.8652.8282.7742.7002.538 

114.8443.9823.5873.3573.2043.0953.0122.9482.8962.8542.7882.7392.7012.6462.5702.404 

124.7473.8853.4903.2593.1062.9962.9132.8492.7962.7532.6872.6372.5992.5442.4662.296 

134.6673.8063.4113.1793.0252.9152.8322.7672.7142.6712.6042.5542.5152.4592.3802.206 

144.6003.7393.3443.1122.9582.8482.7642.6992.6462.6022.5342.4842.4452.3882.3082.131 

154.5433.6823.2873.0562.9012.7902.7072.6412.5882.5442.4752.4242.3852.3282.2472.066 

164.4943.6343.2393.0072.8522.7412.6572.5912.5382.4942.4252.3732.3332.2762.1942.010 

174.4513.5923.1972.9652.8102.6992.6142.5482.4942.4502.3812.3292.2892.2302.1481.960 

184.4143.5553.1602.9282.7732.6612.5772.5102.4562.4122.3422.2902.2502.1912.1071.917 

194.3813.5223.1272.8952.7402.6282.5442.4772.4232.3782.3082.2562.2152.1552.0711.878 

204.3513.4933.0982.8662.7112.5992.5142.4472.3932.3482.2782.2252.1842.1242.0391.843 

214.3253.4673.0722.8402.6852.5732.4882.4202.3662.3212.2502.1972.1562.0962.0101.812 

224.3013.4433.0492.8172.6612.5492.4642.3972.3422.2972.2262.1732.1312.0711.9841.783 

234.2793.4223.0282.7962.6402.5282.4422.3752.3202.2752.2042.1502.1092.0481.9611.757 

244.2603.4033.0092.7762.6212.5082.4232.3552.3002.2552.1832.1302.0882.0271.9391.733 

254.2423.3852.9912.7592.6032.4902.4052.3372.2822.2362.1652.1112.0692.0071.9191.711 

264.2253.3692.9752.7432.5872.4742.3882.3212.2652.2202.1482.0942.0521.9901.9011.691 

274.2103.3542.9602.7282.5722.4592.3732.3052.2502.2042.1322.0782.0361.9741.8841.672 

284.1963.3402.9472.7142.5582.4452.3592.2912.2362.1902.1182.0642.0211.9591.8691.654294.1833.3282.9342.7012.5452. 

077: .0031. 9481. 9041. 8391. 7441. 509504. 0343. 1832. 7902. 5572. 4002. 2862. 1992. 1302. 0732. 0261. 9521. 8951. 8501. 7841. 6871 



Tab.7.7-Tabledesquantilesd’ordre0.99d'unevariabledeFisher'anL etn 2 degesdeliberte 


m=12345 6 7 8 9 10 12 14 16 20 30 oo 

02=1405250005403562557645859 5928 5981 6022 6056 6106 6143 6170 6209 6261 6366 

298.5099.0099.1799.2599.3099.3399.3699.3799.3999.4099.4299.4399.4499.4599.4799.50 
334.1230.8229.4628.7128.2427.9127.6727.4927.3527.2327.0526.9226.8326.6926.5126.13 
421.2018.0016.6915.9815.5215.2114.9814.8014.6614.5514.3714.2514.1514.0213.8413.46 
516.2613.2712.0611.3910.9710.6710.4610.2910.1610.059.8889.7709.6809.5539.3799.020 
613.7510.939.7809.1488.7468.4668.2608.1027.9767.8747.7187.6057.5197.3967.2296.880 
712.259.5478.4517.8477.4607.1916.9936.8406.7196.6206.4696.3596.2756.1555.9925.650 
811.268.6497.5917.0066.6326.3716.1786.0295.9115.8145.6675.5595.4775.3595.1984.859 
910.568.0226.9926.4226.0575.8025.6135.4675.3515.2575.1115.0054.9244.8084.6494.311 

1010.047.5596.5525.9945.6365.3865.2005.0574.9424.8494.7064.6014.5204.4054.2473.909 

119.6467.2066.2175.6685.3165.0694.8864.7444.6324.5394.3974.2934.2134.0993.9413.602 

129.3306.9275.9535.4125.0644.8214.6404.4994.3884.2964.1554.0523.9723.8583.7013.361 

139.0746.7015.7395.2054.8624.6204.4414.3024.1914.1003.9603.8573.7783.6653.5073.165 

148.8626.5155.5645.0354.6954.4564.2784.1404.0303.9393.8003.6983.6193.5053.3483.004 

158.6836.3595.4174.8934.5564.3184.1424.0043.8953.8053.6663.5643.4853.3723.2142.868 

168.5316.2265.2924.7734.4374.2024.0263.8903.7803.6913.5533.4513.3723.2593.1012.753 

178.4006.1125.1854.6694.3364.1023.9273.7913.6823.5933.4553.3533.2753.1623.0032.653 

188.2856.0135.0924.5794.2484.0153.8413.7053.5973.5083.3713.2693.1903.0772.9192.566 

198.1855.9265.0104.5004.1713.9393.7653.6313.5233.4343.2973.1953.1163.0032.8442.489 

208.0965.8494.9384.4314.1033.8713.6993.5643.4573.3683.2313.1303.0512.9382.7782.421 

218.0175.7804.8744.3694.0423.8123.6403.5063.3983.3103.1733.0722.9932.8802.7202.360 

227.9455.7194.8174.3133.9883.7583.5873.4533.3463.2583.1213.0192.9412.8272.6672.305 

237.8815.6644.7654.2643.9393.7103.5393.4063.2993.2113.0742.9732.8942.7812.6202.256 

247.8235.6144.7184.2183.8953.6673.4963.3633.2563.1683.0322.9302.8522.7382.5772.211 

257.7705.5684.6754.1773.8553.6273.4573.3243.2173.1292.9932.8922.8132.6992.5382.169 

267.7215.5264.6374.1403.8183.5913.4213.2883.1823.0942.9582.8572.7782.6642.5032.131 

277.6775.4884.6014.1063.7853.5583.3883.2563.1493.0622.9262.8242.7462.6322.4702.097 

287.5365.4534.5684.0743.7543.5283.3583.2263.1203.0322.8962.7952.7162.6022.4402.064297.5985.4204.5384.0453.7253. 

801/ .6652.5632.4842.3692.2031.805507.1715.0574.1993.7203.4083.1863.0202.8902.7852.6982.5622.4612.3822.2652.098] 
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